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ABSTRACT 


(Distribution  Limitation  Statement  B) 

The  theoretical  problem  of  pulsed  chemical  laser  ignition  by 
very  rapid  infrared  laser  absorption  heating  is  investigated.  It  is  found 
that  for  single-phase  systems,  the  only  hope  for  ever  achieving  energy 
gain  appears  to  lie  in  the  utilization  of  long  or  fast-branching  chemical 
chains  in  the  pumping  reactions.  For  the  hydrogen-halide  system,  the 
usable  chain  length  is  severely  limited  by  the  very  rapid  collisional  de¬ 
excitation  rate  of  the  vibrational  states  by  the  reaction  product  molecules, 
making  positive  energy  gain  very  unlikely.  In  two-phase  systems,  how¬ 
ever,  one  may  try  to  utilize  the  large  amount  of  combustion  energy 
release  from  burning  of  finely-dispersed,  submicron  size  solid  particles 
for  generation  of  additional  atoms  required  in  the  pumping  of  the  laser 
reactions.  The  reduced  input  energy  requirement  for  generating  a  fixed 
number  of  active  atoms  in  such  two-phase  system  may  improve  the 
prospect  for  achieving  positive  energy  gain  in  multistage  pulsed  chemical 


laser  applications. 
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SECTION  I 


summ:  ky 

A  summary  of  the  essential  findings  of  a  6-month  theoretical  study 
on  the  feasibility  of  pulsed  chemical  laser  ignition  by  very  rapid  infrared 
absorption  heating  is  presented.  Specifically,  the  infrared  source  of 
particular  interest  is  the  output  from  another  high-power  laser,  so  that 
the  subject  of  this  investigation  is,  in  effect,  pulsed-chemical  laser 
initiation  by  another  laser.  Detailed  technical  results  from  this  investiga¬ 
tion  are  presented  in  section  II. 

1.  TECHNICAL  PROBLEM  AND  OBJECTIVE  OR  STUDY 

To  make  high-energy  pulsed  chemical  la^er  operations  feasible- 
it  is  necessary  to  obtain  some  multiplication  factor  of  the  input  initiation 
energy.  There  are  currently  several  methods  of  initiating  a  chemical 

laser - flashlamp,  spark  discharge,  and  electron  beam.  All  of  these 

methods  have  potential  limitations  when  applied  to  scaling  of  large 
pulsed  laser  systems.  Among  the  problem  areas  are  efficiency  of  initia¬ 
tion,  volume  uniformity  of  in itiat ion,  speed  of  initiation,  and  the  required 
pulsed  duration  of  the  initiating  system.  One  method  which  has  not  yet 
been  investigated  is  the  initiation  of  a  pulsed  chemical  laser  by  another 
laser.  In  view  of  the  current  rapid  development  of  high-power  chemical 
and  nonchemical  lasers,  and  the  many  desirable  characteristics  of  such 
lasers  (e.  g.  ,  rapid  rise  time,  beam  uniformity,  adaptability  to  multiple 
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staging,  etc.  )  which  suggest  themselves  as  natural  initiation  sources  for 
the  later-stage  power  amplifiers  in  a  pulsed  chemical  laser  system,  the 
feasibility  of  such  initiation  method  merits  exploration  so  that  their 
potential  usefulness  can  be  further  assessed. 

The  objective  of  the  present  investigation  is,  therefore,  to  deter¬ 
mine  the  feasibility  of  initiating  a  pulsed  chemical  laser  by  the  use  of 
another  laser.  The  main  emphasis  of  this  study  is  in  the  identification 
and  definition  of  the  physical  and  chemical  conditions  under  which  such 
initiation  process  would  be  favorable. 

2.  GENERAL  METHOD  OF  APPROACH 

Since  this  is  a  theoretical  study  involving  the  exploration  oi  an 
original  concept,  our  method  of  approach  necessarily  involves  more  than 
what  is  ordinarily  required  in  a  standard  analytical  program.  Thus,  in 
addition  to  literature  search  for  relevant  data  and  related  previous  works, 
it  is  necessary  to  carry  out  conceptual  formulation  involving  consideration 
of  many  alternatives  before  detailed  formulation  and  numerical  computa¬ 
tion  can  begin.  In  areas  where  no  experimental  data  or  prior  knowledge 
exist,  it  is  also  necessary  to  consult  experts  in  related  fields  so  that 
meaningful  extrapolations  from  existing  knowledge  can  be  made  and  the 
probable  uncertainties  in  the  resultant  ordev-of-magnitude  estimates  can 
be  assessed. 
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Because  of  the  extremely  short  duration  of  the  contracted  study, 
opportunity  for  in-depth  investigation  is  severely  limited,  and  breadth 
coverage  of  many  pertinent  subjects  is,  at  best,  highly  selective.  A  best 
effort  has  been  made,  however,  in  striking  a  proper  balance  between 
breadth  and  depth  coverage  within  the  time  and  resources  available  so  that 
the  most  important  features  and  inherent  limitations  of  the  pulsed  chemical 
laser  initiation  problem  can  be  brought  out  and  some  preliminary  conclu¬ 
sions  can  be  drawn  regarding  the  feasibility  of  the  laser  initiation  process. 

3.  TECHNICAL  RESULTS 

From  a  pragmatic  point  of  view,  the  use  of  one  laser  to  initiate  or 
"trigger"  another  laser  will  be  meaningful  only  if  one  or  more  of  the 
following  functions  could  be  accomplished:  (a)  Power  or  energy  amplifica¬ 
tion,  (b)  frequency  shift  to  facilitate  beam  propagation  and/or  energy 
coupling  to  target,  and  (c)  pulse  shaping  for  the  same  purposes.  However, 
because  the  energy  output  per  unit  volume  or  per  mole  of  reactants  from 
current  pulsed  chemical  lasers,  by  whatever  means  of  initiation,  is 
generally  minuscule  in  comparison  with  their  true  potential,  achievement 
of  significant  energy  amplification  have,  therefore,  been  considered  our 
most  important  immediate  objective. 

It  may  be  noted  that  all  chemical  laser  systems  which  have  been 
studied  or  under  active  investigation  as  reported  in  the  published  and 
unpublished  literatures  up  to  the  time  of  commencement  of  the  present 
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research  (May  1971)  were  "single-phase"  systems.  That  is,  the  lasing 
medium  from  which  chemical  energy  was  to  be  extracted  and  converted 
into  coherent  radiation  energy  was  either  aM  in  the  gas-phase  or  all  in  the 
liquid-phase.  In  fact,  the  most  energetic  chemical  laser  systems  dis¬ 
covered  to  date  have  been,  and  apparently  still  are,  purely  single-phase 
systems  involving  only  gaseous  mixtures  (e.  g.  ,  the  hydrogen-halogen 
systems).  For  this  very  reason,  the  proposed  feasibility  study  was  origi¬ 
nally  intended  only  for  single-phase  pulsed  chemical  laser  systems.  This  is 
reflected  in  the  contractual  work  statement  (Section  I).  However,  soon 
after  the  commencement  of  the  contracted  study,  we  discovered  that  artificial 
enhancement  of  infrared  absorption  in  a  laseable  gas  mixture  can,  in  principle, 
be  achieved  by  the  addition  of  finely-disper^pd,  submicron  size  solid 

if 

particles  (or  liquid  droplets)  into  the  gaseous  medium.  Furthermore,  the 
addition  of  such  solid  particles  of  suitably  chosen  chemical  properties  into 
a  laseable  gas  mixture  may  greatly  improve  the  energy  gain  potential  of 
the  resultant  "two-phase  system.  " 

The  consideration  of  more  than  one  phase  of  matter  in  a  lasing 
medium  clearly  represents  an  additional  complication  which  one  wishes  to 
avoid,  if  at  all  possible.  However,  we  feel  that  in  any  search  for  poten¬ 
tially  high  gain  and/or  high-power  laser  systems  in  the  conceptual  stage, 
there  appears  to  be  no  good  a  priori  reason  for  limiting  oneself  to  the 
consideration  of  only  single-phase  systems.  Accordingly,  with  the  con¬ 
sent  of  cognizant  personnel  at  the  Air  Force  Weapons  Laboratory  (AFWL) 

Air  Force  Special  Weapons  Center  (PMRA),  contract  number  F29601-71- 
C-0081,  5  May  1971,  pp.  25-26. 
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and  at  the  Advanced  Research  Projects  Agency/Strategic  Technology  Office 
(AJRPA/STO),  we  have  broadened  the  scope  of  our  theoretical  study  of 
pulsed  chemical  laser  initiation  to  include  the  consideration  of  two-phase 
systems"  --  with  no  increase  in  the  contracted  leve  of  effort 

The  major  conclusions  which  we  have  arrived  at  by  the  end  of  the 
6-month  period  are  as  follows: 

A.  Initiation  of  a  pulsed  chemical  laser  by  another  laser  is 
generally  feasible. 

B.  In  single-phase  pulsed  chemical  laser  systems,  the  only  hope 
for  ever  achieving  energy  gain  appears  to  lie  in  the  utilization 
of  long  or  fast-branching  chemical  chains  in  the  pumping 
reactions. 

C.  For  the  hydrogen-halide  system,  which  is  about  the  only 
chemical  laser  system  that  has  been  demonstrated  to  be  capable 
of  operating  at  high-power  density  in  continuous  wave  (CW) 
application,  the  usable  chain  length  appears  to  be  severely 
limited  by  the  very  rapid  collisional  de-excitation  rate  of  the 
vibrational  states  by  the  reaction  product  molecules  [e.  g.  ,  de¬ 
excitation  of  HF-(v)  by  HF (0 )].  The  chain- branching  rate  in  the 
single-phase  hydrogen-halide  system  also  appears  to  be  too 

* 

slow  to  be  of  any  help. 

D.  If  long  chemical  chains  cannot  be  utilized,  and  the  energy 
required  for  generating  active  atoms  for  the  pumping  of  the 
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lasing  reactions  must  come  from  the  external  energy  source, 
the  laser  output  energy  will  always  be  smaller  than  the  inpuc 
energy  supplied  by  the  external  source  and  hence  a  positive 
energy  gain  will  not  be  possible.  The  fact  that  all  pulsed 
chemical  laser  experiments  using  flashlamp  and  electric  dis¬ 
charge  initiation  in  hydrogen-halogen  mixtures  failed  to  yield 
any  significant  laser  output  energy  in  relation  to  the  initiation 
energy  is,  perhaps,  a  good  indication  of  this  situation. 

E.  In  two-phase  systems,  one  may  try  to  make  use  of  the  large 
amount  of  combustion  energy  release  from  burning  of  finely- 
dispersed,  submicron  size  solid  particles  for  generation  of 
dissociated  atoms  in  the  hot  flame  zone  surrounding  the 
particles.  If  the  combination  of  rate  constants  and  experi¬ 
mental  parameters  were  right,  there  is  a  good  chance  that 
frozen-diffusion  of  dissociated  atoms  from  the  hot  flame  zone 
into  the  bulk  of  the  inter-particle  space  would  lead  to  efficient 
conversion  of  the  dissociation  energy  into  laser  energy  through 
nonchain  rearrangement  reactions. 

F.  The  heat  of  combustion  of  some  metal  and  solid  compound  is 
often  two  orders  of  magnitude  greater  than  the  energy  required 
to  heat  the  same  mass  of  metal  or  compound  up  to  their  respec¬ 
tive  ignition  temperature,  there  exists  a  good  possibility  for 

l 
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achieving  positive  energy  gain  in  the  two-phase  systems. 

G.  Our  present  knowledge  about  the  burning  rate  of  submicron 

size  particles  in  laseable  gas  mixtures  and  about  the  chemical 
kinetic  and/or  de- excitation  rates  of  their  reaction  products 
is  too  incomplete  to  allow  any  quantitative  prediction  oi  two- 
phase  pulsed  chemical  laser  performance.  In  fact,  we  do  not 
even  know  the  initial  chemical  stability  (or  time-dependent 
explosion  limit)  of  any  of  the  energetically  attractive  two- 
phase  mixtures  to  allow  an  assessment  of  their  suitability  for 
premixed  pulsed  laser  operation  at  this  time. 

4.  DOD  IMPLICATIONS 

Existing  DOD/ARPA  pulsed  chemical  laser  programs  at  the 
Aerospace  Corporation,  Avco  Everett  Research  Laboratory,  and  at  other 
contractor  locations  all  involve  single-phase  systems  (i.  e.  ,  all  gas  mix¬ 
tures).  Initiation  mechanisms  for  these  single-phase  systems  are  either 
(a)  flashlamp  initiation,  (b)  electrical  discharge  initiation,  or  (c)  electron 
beam  initiation.  So  far,  no  work  on  two-phase  pulsed  chemical  laser 
appears  to  have  been  done. 

In  view  of  the  positive  energy  gain  potential  indicated  for  the  two- 
phase  pulsed  chemical  laser  system  involving  solid  particle  combustion, 
it  appears  that  a  sufficient  level  of  effort  should  be  mounted  by  DOD  in  the 
exploration  of  this  potential  so  that  its  eventual  development  would  not 
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represent  a  technical  surprise.  The  problem  of  laser  ignition  of  single- 
and  two-phase  pulsed  chemical  laser  systems  also  merits  further  explora¬ 
tion  since  it  may  lead  to  convenient  multiple  staging  of  the  power  (or 
energy)  amplifiers  within  such  systems. 

5.  IMPLICATIONS  FOR  FURTHER  RESEARCH 

The  present  investigation  of  the  feasibility  of  pulsed  chemical  laser 
ignition  by  very  rapid  infrared  absorption  heating  clearly  points  to  the  need 
for  further  research  in  the  following  areas  of  chemical  physics  so  that  the 
energy  gain  potential  of  two-phase  systems  can  be  more  quantitatively 
assessed: 

A.  Chemical  and  thermal  stability  of  two-phase  mixtures  with 
laser  potential. 

B.  Optical  properties  of  solid  particles  which  are  energetically 
attractive  when  used  in  two-phase  chemical  laser  systems. 

4 

C.  Ignition  process  and  the  rate  of  combustion  of  submicron  size 
particles  in  laseable  gas  mixtures. 

D.  Emissivity  of  the  products  of  reaction  in  two-phase  chemical 
laser  systems  involving  solid  particle  combustion. 

E.  Energy  transfer  mechanisms  during  solid  particle  combustion 
in  two-phase  systems. 

F.  Excitation  and  de-excitation  processes  in  heterogeneous, 
nonsteady  chemical  reactions. 
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SECTION  II 

DETAILED  TECHNICAL  RESULT 

1.  INTRODUCTION 

The  potential  of  chemical  lasers  as  efficient,  high  specific  energy 

coherent  radiation  sources  has  been  well  demonstrated  by  the  recent  works 

1  2-5 

of  Airey  and  McKay,  Cool,  Stephens  et  al.  ,  and  Spencer,  Jacobs, 

6-9 

et  al.  Among  the  different  single-phase  chemical  lasers  studied,  the 
deuterium  fluoride  (DF)  and  the  deuterium  fluoride-carbon  dioxide  (DF-CO^ 
systems  within  the  hydrogen-halide  group  appear  most  promising,  both 
from  the  point  of  view  of  energy  conversion  efficiency  and  that  of  atmos¬ 
pheric  propagation.  Of  these  two  systems,  the  pure  DF  laser,  with  strong 
lines  lying  within  the  3.  5  ~  4.  2u  atmospheric  wdndow,  offers  one  possible 
advantage  over  the  DF-CO^  system  in  terms  of  higher  breakdown  threshold 
(by  a  factor  of  about  8  according  to  frequency-square  scaling).  On  the 
other  hand,  the  DF-CO^  system,  being  a  mixed  gas  system  involving 
vibrational  energy  transfer  to  a  chemically  inert  diluent  (i.  e.  ,  CO^)  of 
variable  concentration,  may  offer  the  advantage  of  being  more  flexible  in 
pulse- shaping. 

1-9 

While  all  of  the  recent  works  reported  in  the  literature  have 
been  concentrated  in  the  study  of  CW  chemical  laser  operation  and  per¬ 
formance,  corresponding  effort  in  the  basic  study  of  pulsed  chemical  lasers 
appears  to  be  lagging,  hi  addition  to  the  problem  of  inversion  kinetics. 
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there  is  the  additional  important  problem  of  initiation  in  the  pulsed  opera¬ 
tion,  To  make  high-energy  pulsed  operations  feasible  it  is  necessary  to 
obtain  some  multiplication  factor  of  the  input  energy.  There  are  currently 
several  methods  of  initiating  a  chemical  laser  -  spark  discharge,  flash- 
lamp,  and  electron  beam.  All  of  these  systems  have  potential  limitations 
when  applied  to  scaling  of  iarge  pulsed  laser  systems.  Among  the  problem 
areas  are  efficiency  of  initiation,  volume  uniformity  of  initiation,  speed  of 
initiation  and  the  required  pulsed  duration  of  the  initiating  system.  One 
method  currently  not  being  investigated  is  the  initiation  of  a  pulsed  laser 
with  another  laser. 

The  objective  of  this  study  is,  therefore,  to  determine  the  feasi¬ 
bility  of  initiating  a  pulsed  chemical  laser  by  the  use  of  another  laser. 

The  main  emphasis  of  this  study  is  to  define  the  physical  and  chemical 
conditions  under  which  such  initiation  would  be  favorable. 

The  most  penetrating  theoretical  analysis  of  single-phase  chemical 
laser  operation  applicable  to  pulsed  hydrogen-halide  systems  involving 
chain  chemical  reactions  to  date  appears  to  be  that  of  Oraevskii,  ^  even 
though  the  analysis  was  limited  to  linear  systems  involving  very  few  rate- 
limiting  reactions.  Theoretical  and  experimental  study  of  nonchain  pulsed 
chemical  lasers  in  single-phase,  hydrogen-halide  systems  initiated  by 
flash  photolysis  have  also  been  reported  by  Airey. 
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2.  SOME  FUNDAMENTAL  CONSIDERATIONS  OF  PULSED  CHEMICAL 
LASER  DESIGN  AND  OPERATION 

8  3 

As  explained  by  Spencer  et  aL  and  Cool  et  al.  ,  laser  action  in 
the  continuous  -  flow  DF  and  the  DF-CCL  systems  is  caused  by  population 
inversion  among  the  first  few  vibrational  levels  of  DF,  and  that  between  th 
(10°0)/(02  0)  and  the  (00  1)  levels  of  CO?,  respectively.  For  both  system 
the  chemical  pumping  energy  was  derived  mainly  from  flit:  pair  of  chained 
rearrangement  reactions 


F  +  D  -  DF(v)  1  D 

All  = 

-31.7  kcal /mole 

(1) 

D  +  Y  -  DF(v)  1  F 

A  H  - 

-  1  Ob,  7  Real/ mole 

(2) 

in  which  a-  large  fraction  of  the  respective  chemical  energy  released 
(i.  e,  ,  negative  value  of  the  heat  of  formation  ill)  ended  tip  in  the  vibra¬ 
tional  excitation  of  the  newly  formed  DF  molecules.  In  the  experimental 

£ 

arrangement  of  Spencer  et  al,  ,  only  reaction  (1)  was  involved  and  the 
free  fluorine  atoms  needed  in  the  chained  rearrangement  reaction  was 
supplied  by  thermal  decomposition  of  SF^  ,  using  arc-heated  nitrogen  as 
an  external  heat  source.  In  the  experimental  arrangement  of  Cool  et  al.  , 
reactions  (I)  and  (2)  were  allowed  to  proceed  simultaneously  in  a  mixing 

flow  consisting  of  one  stream  of  partially  dissociated  F  and  another 
stream  of  D^-CO.,  mixture.  Partial  dissociation  of  the  Fn  gas  in  the 


primary  stream  was  accomplished  either  by  a  high-frequency  electrical 
discharge  or  by  an  auxiliary  hypergolic  reaction  ol  the  type 

F  +  NO  -  NO F  4  F  (3) 

lit 

It  is  important  to  note  that  in  the  continuous  flow  chemical  lasers 

2  —  3  ^  (i  —  ^ 

of  the  types  described  by  Cool  et  al,  and  by  Spencer  el  ai,  ,  the 
chemical  pumping  rate  was  actually  limited  by  the  rate  of  molecular  mix¬ 
ing  between  the  active  reactants  when  the  initially  segregated  streams 
were  brought  together,  and  not  by  the  chemical  kinetic  rates  tor  the  cor¬ 
responding  homogeneous  gas-phase  reactions.  The  de- excitat ion  rate, 
on  the  other  hand,  was  controlled  by  vibrational  relaxation  ol  the  respec¬ 
tive  lasing  molecular  species. 

The  use  of  heterogeneous  stream  mixing  as  a  means  lor  controlling 
the  chemical  reaction  and  laser  pumping  rates  by  the  above  cited  investi¬ 
gators  were  probably  necessitated  by  the  fact  that  a  premixed  gas  con¬ 
taining  substantial  concentrations  ol  very  active  chemicals  like  F  and 
has  a  tendency  of  being  easily  detonated  by  a  fast  chemical  chain 
formed  by  a  close  coupling  between  reactions  (1)  and  (2),  and  is  therefore 
too  difficult  to  handle  in  a  continuous  flow  system. 

Because  the  maximum  output  power  per  unit  active  volume 
for  any  laser  system  is  directly  proportional  to  the  population 
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inversion  density,  an  obvious  direction  for  increasing  the  chemical  laser 
output  power  density  is  to  increase  the  reactant  mass  density  within  the 
mixing  reactor  while  trying  to  maintain  a  fixed  inverted  population  ratio 
between  the  upper  and  the  lower  lasing  states.  However,  rapid  mixing  of 
initially  heterogeneous  gas  streams  of  high-mass  densities  generally 
necessitates  the  use  of  high  intensity  turbulent  flows.  Furthermore, 
operation  at  high  gas  densities  also  tends  to  shorten  the  vibrational  relaxa¬ 
tion  time  for  de-inversion  of  the  laser  levels.  This  tends  to  limit  the 
available  streamwise  dimension  of  the  active  region.  Thus,  one  of  the 
major  problems  in  future  development  of  high  power-density,  continuous- 
wave  chemical  lasers  of  good  optical  quality  may  well  lie  in  the  area  of 
controlling  the  spatial  homogeneity  and  temporal  stability  of  high  Reynolds 
number  compressible  turbulent  flows  with  strong  chemical  reactions  --  a 
well  known  unsolved  problem  in  contemporary  gas  dynamics! 

In  the  conceptual  design  of  high-energy  pulsed  chemical  lasers 

-  4 

of  relatively  short  pulse  duration  (say,  of  the  order  ol  10  sec  or  less), 
the  introduction  of  convective  flow  motion  appears  to  offer  no  particular 
advantage,  except  for  the  purpose  of  gas  sample  replacement  within  the 
active  volume  for  repetitive  operations.  The  argument  is  simply  that: 

(a)  if  a  premixed  reactant  mixture  (e.  g.  ,  a  stoichiometric  mixture  of 
and  gases  at  standard  temperature  and  pressure)  were  used  in  a  con¬ 
tinuous-flow  arrangement  through  the  active  laser  volume,  the  danger  of 
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accidentally  detonating  the  entire  reactant  storage  and  supply  system  by 
the  pulsed  chemical  reactions  would  he  too  great;  (b)  if  an  uii-premixod 
flow  system  were  used  to  bring  the  reactants  into  the  active  laser  volume 
in  a  continuous  fashion,  molecular  mixing  of  tin'  reactants  may  not  have 
time  to  reach  completion  during  the  intended  pulse  period,  resulting  in  a 
loss  of  effective  lasing  volume  or  an  undesirable  stretching  out  of  the  laser 
pulse,  or  both.  Furthermore,  one  would  have  to  lace  the  same  problem  of 
controlling  the  homogeneity  of  the  reacting  mixture  as  in  the  case  of  a 
continuous-wave  chemical  laser  considered  in  the  preceding  paragraph; 

(c)  if  the  reactants  were  brought  into  the  active  laser  volume  in  an  inter¬ 
mittent  manner,  then  there  appears  10  be  no  compelling  reason  tor  initiat¬ 
ing  the  pulsed  chemical  reaction  belorc  most  ot  the  residual  fluid  motions 
have  subsided. 

The  above  argument  then  leads  one  to  the  tentative  conclusion  that 
a  high-energy  pulsed  chemical  laser  would  most  likely  consist  of  the 
following  essential  elements: 

(a)  An  explosion-proof  pulsed  chemical  laser  reactor,  with  infra¬ 
red  window  (or  windows)  and  other  necessary  optical  elements 
for  coupling  out  the  laser  energy.  The  internal  surfaces  of  the 
reactor,  of  course  must  be  suitably  passivated  so  as  to  reduce 
the  probability  of  prematurely  detonating  the  reactants  by 
surface-catalysed  reactions  to  a  minimum. 
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(b)  Explosion-proof  and  last-acting  gas  filling  and  evacuation 
systems  fur  loading  of  reactants  into,  and  removal  of  reaction 
products  out  of,  the  pulsed  reactor  within  certain  desired  pulse 
repetition  period.  If  the  reactants  were  not  to  be  premixed, 
the  gas  filling  system  must  be  suitably  designed  to  ensure 
uniform  mixing  of  the  reactants  within  the  allocated  time 
period  (i.  e.  ,  some  fraction  of  live  pulse  repetition  period). 

(c)  Provisions  lor  igniting  the  pulsed  chemical  reactions 
uniformly  over  the  active  volume  of  the  laser  reactor. 

The  provision  of  items  (a;  and  (b)  appears  relatively  straightforward. 
It  is  only  a  matter  of  degree  to  which  various  mechanical  per forinance 
parameters  demanded  of  the  laser  system  can  actually  be  satisfied  (e.  g.  , 
maximum  operating  pressure,  maximum  pulse  repetition  rate,  etc.  ).  The 
last  provision  (c),  on  the  other  hand,  is  expected  to  be  the  most  crucial. 

It  is  precisely  the  suggestion  of  a  new'  method  for  homogeneous  ignition  of 
pulsed  chemical  lasers,  together  with  an  exploration  ul  its  potential,  that 
provided  the  initial  motivation  for  the  present  study. 

3.  METHODS  OF  INITIATION 

Eet  us  now  consider  the  crucial  question  of  uniform  ignition  of  a 
detonable  mixture  like  or  H  -  E?  ,  in  a  pulsed  chemical  laser  sys¬ 

tem.  The  conventional  mothud  of  point-ignition  to,  g.  ,  use  of  spark-plugs, 
as  in  the  case  of  internal  combustion  engines)  is  bound  to  be  unsatisfactory. 
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}  rom  (’li'iiiontary  physics  considerations,  there  appear  to  be  only  two 
general  methods  for  effecting  simultaneous  (or  "nearly"  simultaneous 
when  measured  against  the  characteristic  time  scale  of  the  chemical 
reaction)  ignition  over  an  extended  volume,  namely:  (a)  pulsed  electro¬ 
magnetic  wave  absorption;  and  (b)  fast  particle  injection.  We  will  briefly 
consider  the  essential  difference  between  these  two  methods  here: 
a.  Initiation  by  Pulsed  Idled romagnetic  Wave  Absorption 

This  method  can,  in  principle,  be  used  whereve  .here  exists  a 
strong  absorption  coefficient  in  the  laseable  mixture  over  the  frequency 
range  of  (he  source  of  electromagnetic  radiation.  The  resultant  spatial 
homogeneity  of  the  ignited  mixture  is  expected  to  depend  on  the  initial 
spatial  homogeneity  of  the  radiation  source,  and  of  the  mixture  absorption 
coefficient.  Obviously,  in  order  to  achieve  nearly  homogeneous  volume 
absorption  from  an  external  radiation  source  (or  sources)  introduced 
through  the  boundary  surfaces  of  the  reactive  \  olume,  the  absorption 
coefficient  must  not  be  too  strong  so  as  to  cause  excessive  attenuation  of 
the  incident  radiation.  On  the  other  hand,  in  order  10  achieve  simultaneity 
of  ignition  over  the  reaction  volume,  the  combination  of  absorption  coeffi¬ 
cient  and  incident  radiation  intensity  must  be  sufficiently  high  so  that  the 
characteristic  time  for  mixture  heating  remains  short  in  comparison  with 
the  characteristic  lime  for  the  laser-pumping  chemical  reaction. 


According  to  the  nature  of  the  source  of  electromagnetic  radiation, 
this  method  can  further  be  subdivided  into  two  or  more  classes; 

(1)  Flashlamp  Ignition.  The  most  common  meLhod  for  initiating 
a  pulsed  chemical  laser  generally  employs  an  electronic  flashlamp  as  the 
source  of  electromagnetic  radiation.  Since  the  emission  spectrum  of  high 
intensity  flashlamps  usually  peaks  in  the  near  ultraviolet,  the  incident 
radiation  can  be  used  for  direct  generation  of  chain  centers  (e.  g.  ,  halogen 
atoms  in  hydrogen-halide  systems)  through  phutodissociation  of  the  reac¬ 
tant  molecules  in  most  of  the  chemical  laser  systems  studied  to  date. 

This  method  of  initiation  is,  therefore,  also  called  the  method  of  flash 
photolysis  initiation. 

Since  the  process  of  flash  photolysis  is  capable  of  selective  dis¬ 
sociation  of  the  desired  chemical  species  through  proper  matching  of  the 
flashlamp  emission  spectrum  and  the  species  photodissociatiun  spectrum, 
the  principal  advantage  of  this  method  ol  initiation  is  in  its  ability  to  gen¬ 
erate  a  relatively  high  cone entrat ion  of  active  atoms  without  excessively 
heating  up  the  laseable  mixture.  A  secondary  advantage  of  this  method  is 
that  electronic  flashlamps  are  readily  available  and  ol  relatively  low  cost 
(say,  in  comparison  with  a  laser  light  source  of  comparable  pulse  energy). 

The  principal  disadvantage  of  this  method  of  initiation,  on  the  other 
hand,  is  the  limitation  of  radiation  flux  density  (or  light  intensity)  that  can 
be  obtained  from  any  electronic  flashlamp  operating  at  a  finite  plasma 
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temperature.  This  limitation  imposes  an  upper  bound  on  the  rate  of 
photolysis,  and  also  a  requirement  for  large  window  areas  for  coupling  in 
of  the  initiation  energy  to  the  reactive  volume.  A  secondary  disadvantage 
of  this  method  of  initiation  is  that  flashlamps  are  inherently  inefficient 
devices  for  conversion  of  electrical  input  energy  into  light  output  energy 
(say,  in  comparison  with  the  10%  or  higher  conversion  efficiency  that  has 
been  demonstrated  for  some  infrared  electrical  and  chemical  lasers). 
Such  low  conversion  efficiency  tends  to  increase  the  size  and  mass  of  the 
electrical  energy  storage  (e.  g.  ,  capacitor  bank  or  flywheel  generator) 
required  to  deliver  a  given  initiation  pulse  energy  --  which  is,  of  course, 
a  serious  disadvantage  for  any  portables  system! 

(2)  Laser  Ignition.  This  is  the  subject  of  the  present  investiga¬ 
tion  about  which  we  have  much  more  to  say  later.  At  this  point,  it  suf¬ 
fices  to  say  that,  potentially,  laser  ignition  offers  the  same  principal 
advantage  as  does  flashlamp  ignition,  but  without  many  of  the  latter's 
inherent  limitations.  Furthermore,  since  the  input  and  output  radiation 
arc  of  the  same  nature  (i.e.,  both  are  coherent  laser  radiation),  this 
method  of  initiation  readily  renders  itself  to  the  possibility  of  multiple 
staging  in  which  the  laser  output  from  one  power  amplifier  stage  can  be 
used  both  as  the  initiation  energy  source  and  as  the  input  coherent  beam 
for  further  amplification  in  the  next  power  amplifier  stage. 
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b.  Initiation  by  Electronic  Processes 

Electron- molecule  collisions  at  sufficiently  high-electron  energy 
generally  lead  to  dissociation  as  well  as  vibrational  and  electronic  excita¬ 
tion  of  the  molecules.  Thus,  passage  of  an  electric  current  is  a  well  proven 
method  for  initiating  a  chain  chemical  reaction.  However,  achievement  of 
spatial  homogeneity  for  the  resultant  chemical  reaction  is  not  an  easy 
matter  since  the  local  concentration  of  reaction  chain  centers  (e.  g.  ,  dis¬ 
sociated  atoms)  would  be  proportional  to  the  local  electron  current  density, 
and  the  achievement  of  uniform  current  density  distribution  in  a  gaseous 
medium  is  not  at  all  trivial,  especially  at  high  gas  pressures. 

According  to  the  way  in  which  the  initiating  electron  current  is 
introduced  into  the  reaction  volumne,  this  method  can  further  be  subdivided 
into  two  or  more  classes: 

(1)  Pulsed  Electric  Discharge  Ignition.  In  recent  development  of 

12-17 

transversely  excited  atmospheric  (TEA)  CO^  ^asorsi  it  lias  been 

demonstrated  that  relatively  uniform  pulsed  electrical  discharges  in  diatomic 
and  polyatomic  gas  mixtures  at  pressures  comparable  to,  or  somewhat 
higher  than,  the  standard  atmospheric  pressure  (760  torr)  can  lie  achieved 
either  through  the  use  of  multiple,  individually  ballasted  electrodes,  or 
through  the  use  of  double-discharge  with  carefully  contoured  electrodes. 

While  the  degree  of  spatial  homogeneity  so  achieved  lias  been  found  to  be 


satisfactory  for  electrical  pumping  of  chemieally-passive  molecular  lasers, 
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one  does  not  know  at  this  time  it'  it  would  also  be  adequate  for  the  initia¬ 
tion  of  chemical  lasers.  At  any  rate,  for  the  initiation  of  those  chemical 
lasers  involving  chain  chemical  reactions,  the  discharge  current  must  be 
uniformly  distributed  over  the  bulk  of  the  reaction  volume  in  a  time  scale 
that  is  short  in  comparison  with  the  characteristic  chain  propagation  time 
so  as  to  avoid  the  formation  of  detonation  waves.  This  tends  to  impose  a 

more  stringent  requirement  on  the  discharge  geometry  and  method.  (For 

lh 

example,  helical  discharges  which  leaves  a  major  fraction  of  the  reac¬ 
tion  volume  devoid  of  discharge  current  a.re  bound  to  be  unsatisfactory. 

1  7 

The  use  of  double  discharge  may  also  be  inhibited  by  prematured 
detonation  considerations.  ) 

(2)  Electron  Beam  Ignition.  The  possible  vise  of  high  - energy 
electron  beam  as  a  means  for  uniform  pumping  of  CO  -N  lasers  at 
moderately  high  pressure  was  suggested  by  the  principal  investigator  to 
one  of  his  graduate  students  in  1966,  The  advantage  of  using  high- 
energy  electrons  is  that  multiple  scattering  of  the  primary  beam  and  rapid 
spreading  of  the  secondary  electron-ion  pairs  in  the  target  gas  tend  to 
make  the  current  distribution  much  more  uniform  than  would  be  possible 
in  ordinary  electrical  discharges.  The  only  disadvantage  of  using  high- 
energy  electron  beam  is  that  relatively  high  voltage  would  be  needed  in 
the  generation  of  such  beam,  and  that  relatively  thin  windows  made  up  of 
low-Z  (i.c,  ,  low  atomic  number)  material  would  have  to  be  employed  as 
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partition  between  the  beam  generation  chamber  and  the  laser  gas  excita¬ 


tion  chamber. 

The  predicted  superiority  of  high  energy  electron  beam  over  that 
of  conventional  electrical  discharge  as  a  means  for  the  pumping  of  electri¬ 
cal  lasers  ha. 5  since  been  demonstrated  by  recent  works  at  the  AVCO 
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Everett  Research  Laboratory  and  elsewhere.  Purely  from  the  point 
of  view  of  spatial  homogeneity,  one  may  expect  that  the  use  of  high-energy 
electron  beam  woaild  also  be  superior  to  the  use  of  conventional  electrical 
discharge  for  the  purpose  of  igniting  a  chemical  laser.  However,  on 
absolute  term  and  in  actual  application,  there  are  other  factors  to  be 
considered.  For  example,  it  is  not  clear  at  this  point  to  what  extent 
would  the  presence  of  free  electrons  be  detrimental  to  the  maintenance  of 
vibrational  population  inversion  in  an  ignited  chemical  laser  mixture. 
(According  to  the  principle  of  detailed  balance,  electron  impact  processes 
that  are  efficient  in  vibrational  excitation  would  also  be  efficient  in  vibra¬ 
tional  de  •excitation.  )  Such  questions,  of  course,  cannot  be  answered 
without  further  investigation  --  and  investigation  into  such  problems  is 
quite  outside  the  scope  of  the  present  study. 

4.  LASER  1C  Nil  ION  OF  SINGLE- PHASE  SYSTEMS 

Let  us  now  turn  to  the  main  subject  of  the  present  investigation, 
namely,  laser  ignition.  In  this  subsection,  we  will  consider  only  single- 
phase  chemical  laser  systems.  La  particular,  we  will  focus  our 
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attention  on  the  hydrogen-halide  systems  since  such  systems  appear  to  be 
the  only  ones  that  have  been  demonstrated  to  be  capable  of  yielding  rela¬ 
tively  high  power  density  among  all  presently  known  chemical  lasers, 

a.  General  Nature  of  Interaction  between  an  Intense  Laser  Beam 
and  an  Absorbing  Gas 

Interactions  between  a  laser  beam  and  an  absorbing  medium  are 
governed  by  the  electromagnetic  wave  propagation  equations;  the  mass, 
momentum,  and  energy  conservation  equations  of  the  medium;  the  equation 
of  state  which  relates  the  local  stress  (including  the  hydrostatic  pressure) 
to  the  local  strain  and  other  the  nondynamic  state  variables  (such  as  mass 
density  and  translational  tempo  tutu  re  )  oi  the  medium;  the  constitutive 
equations  which  relate  tnr  local  macroscopic  ole<  1  rnnugi.i  tie  properties  of 
the  medium  (such  as  the  refractive  index,  the  absorption  coefl  icients,  etc.  ) 
to  the  local  thermodynamic  state  of  the  medium;  the  relaxation  equations 
(including  all  relevant  chemical  kinetic  rate  equations);  and  also  the 
boundary  conditions  for  the  interaction  medium.  The  general  problem  is 
therefore  very  complex,  and,  in  fact,  open-ended. 

For  the  purpose  of  studying  the  chemical  laser  initiation  problem  in 
a  gaseous  medium,  one  may  ignore  all  quantized  field  effects  and  adopt  the 
somiclassical  treatment  of  electromagnetic  radiation  as  an  approximation. 
Furthermore,  the  processes  of  absorption  and  refraction  can  be  considered 
independently,  with  the  ray  curvature  determined  entirely  by  refractive 


index  gradient  as  in  a  lossless  dielectric  medium.  The  refractive  index 


gradient,  in  turn,  is  determined  by  the  non-uniform  heating  effects  caused 
by  intial  beam  inhomogeneity,  and  non-uniformity  in  the  spatial  distribution 
of  absorption  coefficient,  if  any.  This  is  reasonable  since  absorption  occurs 
along  the  path  length  and  has  little  effect  on  bending  the  rays  except  indirectly 
through  localized  heating,  which  changes  the  refractive  index  distribution. 

(1)  The  Governing  Equations.  Within  the  above  stated  approxima¬ 
tion,  the  radiation  field  associated  with  the  laser  beam  can  be  represented 
by  a  slightly-curved,  nearly-monochromatic  electromagnetic  wave  of  local 
electric  and  magnetic  field  vectors 


E  =  E  e 
o 


-  iu)  t 


(4) 


H  =  H  e 
o 


-in)  t 


(5) 


The  vector  field  amplitudes  E  ,  II  ,  in  turn,  are  governed  by  the  source- 

o  o 


free  Maxwell's  equations,  such  that 
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vxH  +  ik  e  E  =  0 

O  o  o 


(6) 


V  x  E  -  ik  u  II  =  0 
o  o  m  o 


V  .  u  11=0 
m  o 


v  •  £  E  -  0 
o 


(7) 

(«) 

(9) 


where  e  and  u  are  the  local  dielectric  constant  and  magnetic  perrnis- 
m 

sivity  as  defined  by  the  constitutive  equations 
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-  electric  displacement 


(10) 


D  =  c  E 

o 


Bo  =  MmHD  =  magnetic  induction 


and  k  =  2 tt /\  --jj/c  is  the 


magnitude  of  the  vacuum  propagation  vector. 


After  substitution  of  Eq.  (7)  into  (6)  and  applying  (9),  one  gets 


V  V-m  kn  +  V  lo8  U  m  Xpx  K  ) 


+  V  (E^  •  v  log  e)  =  0  (12) 


l'or  propagation  through  gaseous  media,  it  is  reasonable  to  assume 
u  m  =  constant  =  1>  such  that  (12)  becomes 

r/  2  j*  2  2  -*  -♦  -* 

V  Eo+n  ko  EC  4  27<VV  ">  =  0  (13 

where  n  =  Vc  is  the  ordinary  index  of  refraction. 

In  a  gas,  n  =  1  +  i\  n  where  An  «  1,  and  hence  the  ratio  of  the 
third  term  to  the  second  term  in  Eq.  (13)  is  of  the  order  of  magnitude2 3 

E  k  An/f 

O  o  7. 

—  2  A  i]  /  In  1<  D  \ 

En"  k"  . . o"  <14' 

o  O 

where  l  is  the  smallest  scale  of  variation  of  the  refractive  index.  It  can 
always  be  assumed  that  \  «  I  ,  hence  Zhn/n^k^!  «  l  and  therefore 
Eq.  (13)  reduces  to 


v  E  +  n  k 2  E 
o  o  o 


Next,  let  us  assume  that 
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E  (r,  z)  =  A(r,  z)  e 
o 


iTk  z  -  u)  t  +  k  P(r,  z)J 
L  o  o  J 


:i6i 


is  a  component  of  the  laser  beam  electric  field  where  A(r,  z)  and  p  (r,  z) 

are  slowly  varying  functions  of  the  axial  coordinate  z  and  rapidly  varying 

functions  of  the  radial  coordinate  r.  Substituting  Eq.  (16)  into  (15), 

equating  real  and  imaginary  parts,  and  making  the  eikonal  approximation 

with  k  -*  oo  one  obtains 
o 


(17) 


3  A2 

aA2 

.  2 

/  q  ,  3  q 

— —  -{■ 

q  - — 

+  A 

-  +  ^ 

3  z 

o  r 

\  r  3  r 

3  q  , 

3q 

1 

3  n  3  T 

,~z  + 

q  r —  = 

o  r 

n 

3  T  3  r 

(18) 


dP 

where  q  =  —  .  The  condition  q  =  0  imolies  that  phase  fronts  are  moving 
dr 

in  the  z  direction  or  that  the  ray  paths  are  not  diverging.  Equations  (17) 
and  (18)  show  that  the  amplitude  and  phase  of  the  laser  radiation  are  not 

independent.  The  effect  of  absorption  can  be  taken  into  account  by 
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modifying  Eq.  (17)  to  include  the  absorption  term  a  A  ,  such  that. 


hA2  3  A2  2 

+  q  - -  +  A 


5  z 


'3  r 


(? ♦  B) 


h  a  A  -  o 


(19) 


where  a  is  the  absorption  coefficient  which  may  be  a  function  of  the  spatial 
coordinates.  Thus,  for  the  simplest  case  of  a  nondivergent  beam  with  q  =  0, 
Eq.  (19)  reduces  to 
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0 


(20) 


+ 


which  is  the  radiative  transfer  equation  applied  to  a  parallel  beam 
propagating  in  the  z  direction  with  no  source  term. 

Neglecting  convection  and  relaxation  effects,  the  local  gas  tempera¬ 
ture  T  is  related  to  the  electric  field  amplitude  by  the  energy  equation  in 
a  constant  volume  absorption  process 


z 

where  I  =  I(r,  z)  =  local  beam  intensity  =  (c/8nKe(ExlT))  -  c/8n  A  with 

the  {  }  brackets  indicating  a  time  average,  C  is  the  specific  heat  at 

constant  volume,  and  K  is  the  heat  conduction  cuetlieient.  ot  the  gaseous 

medium.  Thus,  if  the  absorption  coefficient  a  were  known,  neglecting 

convection  and  relaxation  effects,  Eqs.  (18),  (19).  and  (21)  constitute  three 

2 

equations  and  three  unknowns  q,  A  ,  and  T.  Those  equations  are  difficult 

21.  22 

to  solve  '  in  general  and  therefore  the  usual  procedure  is  to  solve  the 
equations  assuming  q  =  0  (i.  e.  ,  negligible  beam  divergence)  and  then 
determine  thermal  defocusing  effects  in  an  iterative  manner.  This  is  a 
reasonable  procedure  only  when  the  beam  divergence  effect  is  weak  (i.  e.  , 
either  weak  absorption,  or  small  initial  beam  inhomogeneity,  or  both). 
When  the  beam  divergence  effect  is  strong,  Eqs.  (18),  (19),  and  (21)  must 
be  solved  simultaneously  together  with  the  appropriate  boundary  conditions, 
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including  wall  reflection  and  scattering  to  the  extent  that  they  are 
present. 

{?,)  Thermal  Defocusmg.  Once  the  temperature  distribution  is 
determined  as  indicated  above  for  q  -  0,  the  ray  paths  of  the  laser  beam 
can  be  calculated  within  the  eikonal  approximation  assuming  propagation 
through  a  lossless  dielectric  medium  with  a  known  temperature  dependence 
of  the  refractive  index  n=n(T).  The  problem  would,  of  course,  be  much 
more  complicated  if  relaxation  effects  were  present.  In  the  case,  the 
refractive  index  would  also  be  a  function  of  time,  and  so  would  be  the 
absorption  coefficient  a. 
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The  eikonal  approximation  to  Maxwell's  equations  assumes  that 
the  time  independent  part  of  the  electric  ana  magnetic  field  is  described  by 

ik  Q  (r) 

E  -  £(r)  e  °'  (22) 

o 

ik  C  (r) 

H  -  K(r)  e  °  (2  3) 

o 

where  Q  (r)  is  the  scalar  eikonal  term. 

If  Eqs.  (22)  and  (23)  are  substituted  into  Eqs.  (u),  (7),  and  (8)  and 
going  to  the  geometric  optics  limits  where  k  -*  a.,  one  obtains  the  eikonal 
equation. 

(VC)2  =  n2(r)  (24) 
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Also,  S  -  time  averaged  Poynting  vector  -  c/8  ti  Re(£  xK)  =  c/Su  (G  •  fi  ' )  7  C  . 
Thus,  in  the  geometric  optics  limit,  the  gradient  of  the  eikonal  has  the 
same  direction  as  the  Poynting  vector  or  the  energy  flux  vector.  Therefore, 


(25) 


where  §  is  a  unit  vector  in  the  direction  of  S  and  r(s)  is  a  positive  vector 
along  the  ray  path.  As  the  laser  beam  propagates  through  an  absorbing  gas, 
self  heating  effects  create  radial  temperature  gradients.  These  temperature 
gradients  induce  radial  gradients  in  the  refractive  index  which  cause  the 
beam  to  diverge  and  reduce  the  centerline  intensity.  If  wall  reflection  and 
scattering  were  absent,  the  reduction  in  intensity  due  to  divergence  of  the 
rays  is  determined  from  the  condition 

V  •  S  =  V  •  (if)  =  V  •  f(IV(;/n)J  =  0  (26) 

where  I  is  the  local  intensity.  The  intensity  ratio  at  any  two  points  of  a 
diverging  ray  is  therefore  given  by 


h'h  ‘  'W  esp 


(27) 


The  local  curvature  of  a  ray  is  related  to  the  variation  in  the  refractive 
index  by  use  of  Eq.  (25).  Thus, 
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the  curvature 


where  0  is  the  angular  deflection  of  the  ray.  For  the  special  case  where 
the  radial  gradients  in  the  refractive  index  are  much  greater  than  in  the  z 
direction,  Eq.  (.31)  reduces  to 


1  hi 
-  _ —  d  /. 

n  c  r 


(32) 


gives  the  angular  divergence  of  a  ray  propagating  in  the  z  direction.  To 

first  order,  therefore,  Eqs.  (27)  and  (32)  determine  thermal  defocus  in g 
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effects  as  discussed  by  Gordon  et  al.  and  D.  C.  Smith  for  a  given 
temperature  distribution. 

An  alternate  formulation  which  eliminates  the  iteration  procedure 
is  to  use  the  energy  equation  for  the  electromagnetic  field.  Tints,  the 


energy  equation  is 


v  •  (IS)  4  a  I  = 


+  a  I  =  0 


(33) 


+  lv 


which  reduces  to  the  radiative  transfer  equation 


dl 

ds 


4 


a  I  -  0 


(34) 


for  the  case  S=  constant.  The  term  V  *5  means  that  the  intensity  is 
reduced  by  beam  divergence  as  well  as  absorption.  The  divergence  effect 
is  determined  from  the  eikonal  approximation 


v'  n 


(35) 


and  the  energy  equation  for  the  fluid  such  as  Eq.  (21)  when  convection  effects 
are  neglected.  The  refractive  index  n-n(T)  is  assumed  a  known  function. 
Then  Eqs,  (21),  (33), and  (35)  give  five  equations  and  live  unknowns  1,  I 
and  §  which  must  be  solved  simultaneously.  This  is  more  difficult  than 
the  previously  discussed  iteration  procedure  where  there  are  only  two 
unknowns  I  and  T. 

Using  the  above  formulation,  numerical  calculations  on  the  time- 

dependent  temperature  and  beam  intensity  distributions  in  the  strong  thermal 

interaction  region  of  a  suddenly  turned  on  CO^  laser  beam  propagating 

through  a  pure  CO^  gas  confined  in  a  cylindrical  chamber  of  inside  diameter 

much  larger  than  Lhat  of  the  laser  beam  with  a  Gaussian  initial  intensity 

18  Z  8 

distribution  have  beta  carried  out  by  R.  A.  Chodzko,  ’  using  the  finite- 
dj  Herein  o  method.  In  Chodzko's  calculations,  relaxation  effects  were 
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neglected,  and  the  Imnperature-dependent  absorption  coefficient  was  cal¬ 
culated  according  to  the  assutnption  of  local  equilibrium  population  of  the 
CO^  vibrational  and  rotational  states.  While  not  directly  applicable  to  the 
present  problem  of  interest  (where  the  reactive  gas  container  diameter  is 
likely  to  be  comparable  to,  or  even  somewhat  smaller  than,  the  incident 
beam 'diameter  in  order  to  avoid  excessive  nonuniform  heating),  Chodxko's 
results  clearly  showed  that  at  high  beam  intensities,  nonlinear  healing 
effects  caused  by  the  temperature-dependence  of  the  absorption  coefficient 
are  often  dominant  in  determining  the  time-dependent  temperature  and  beam 
intensity  distributions  wit.iin  the  absorption  region. 

1  8 

The  finite-difference  code  developed  by  Chodzko  can  readily  be 

adopted  for  the  calculation  of  transient  heating  of  any  chemically  passive 

gas  mixture  of  known  a.  (T)  and  n(T)  for  any  desired  container  diameter 

and  initial  beam  intensity  distribution.  With  suitable  modifications,  it  can 

also  be  used  for  the  calculation  of  other  interaction  effects,  such  as 

Z  9  3  0 

optical-acoustic  coupling  ’  and  internal  pressure  wave  excitation. 
However,  these  have  not  been  carried  out  in  the  present  study  due  to  the 
very  limited  time  and  resources  available. 

b.  Gas-phase  Absorption  Processes  Resulting  in  Dissociation 

When  a  laser  beam  is  used  as  an  energy  source  for  the  initiation  of 
a  chemical  laser,  the  must  ideal  absorption  process  would  be  one  in  which 
most  of  the  absorbed  energy  is  spent  in  the  generation  of  active  atoms 


which  contribute  to  the  chemical  pumping  of  the  laser  states.  This  means 


not  only  that  very  little  energy  would  be  wasted  in  translational/ rotational 
heating  and  in  other  unwanted  internal  excitations,  but  also  that  most 
undesirable  interaction  effects  such  as  thermal  defocusing,  optical -acoustic 
coupling,  nonuniform  heating,  etc.,  would  accordingly  be  minimized. 

In  gas-phase  absorption,  there  are  two  known  processes  which  could 
approach  the  above  cited  ideal  when  the  incident  laser  beam  frequency  is  in 
the  right  range. 

(1)  Photodissociation.  It  is  well  known  that  most  halogen  mole¬ 
cules,  such  as  1  Cl^  ,  ^2'  ’  *iavc  strong  dissociation  continua  in  the 
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near  ultraviolet.  Thus,  high -power  pulsed  lasers  ra<  log  in  the  near 
ultraviolet  region  could  conceivably  be  utilized  as  efficient  initiation  sou rcea 
for  hydrogen-halide  lasers  through  direct  pholodissoc  lalion  uf  the  halogen 
molecules.  However,  recently  discovered  coherent  quantum  effects,  such 
as  self- induced  transparency  (SIT)?  ’  may  limit  the  applicability  oi 
photouis sue iatiou  of  such  gases  by  coherent  laser  sources  to  relatively 
high  gas  pressures  or  relatively  low  incident  beam  intensity. 

Even  if  coherent  quantum  effects  does  nut  turn  out  to  be  a  serious 
impediment  lu  the  piloted issocialion  process,  the  feasibility  ui  laser  igni¬ 
tion  through  direct  photodissociation  would  still  hinge  on  the  availability  of 
efficient,  high-energy  pulsed  laser  sources  operating  in  the  near  ultraviolet 
part  of  the  spectrum,  as  well  as  on  the  kinetics  problem  common  to  all 
methods  of  initiation. 
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(2)  Cascade  Excitation  of  Vibrational  Modes.  In  contrast  to  the 
uncertain  availability  of  efficient,  high-energy  pulsed  near  ultraviolet 

lasers,  recent  progress  in  gas  laser  technology  has  made  near-term  avail- 
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ability  of  efficient,  high-energy  pulsed  infrared  lasers  a  certainty.  Thus, 

in  the  laser  ignition  problem,  it  would  lie  most  fruitful  to  look  for  selective 

infrared  absorption  processes  which  may  lead  to  preferential  dissociation 

of  the  participating  molecules.  The  cascade  excitation  dissociation  of  BCl^ 
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molecules  by  CO,,  radiation  as  recently  observed  by  Karlov  et  al.  appears 
to  be  just  one  of  such  processes.  Even  though  it  is  not.  clear  at  Lhis  time 
if  such  selective  dissociation  process  can  proceed  in  laseable  gas  mixtures, 
its  possible  usefulness  as  an  initiation  step  certainly  deserves  further 
exploration. 

c.  Prospect  of  Pulse  Energy  Gain  from  Nonchain  Reactions 

Pulsed  chemical  laser  actions  involving  nonchain  reactions  of  the 

type 

Cl  4  HBr  -  MCI  H  Dr  H  15.5  kcal/mole  (36) 

Cl  +  HI  -  HC1  H  I  H  32.  0  kcal/mole  (37) 


11  35,36 

have  been  reported  by  Airey  and  by  Polanyi  et  al.  In  these  experiments, 


Tectronic  flashlamps  were  used  as  the  energy  source  for  generation  of  the 


active  atoms  through  photodissociation  of  Cl^.  In  such  lionchain  reactions, 
the  laser  output  energy  that  can  lie  extracted  from  each  mule  of  Cl  atoms 


consumed  is  at  most  some  fraction  of  the  exothermic  it y ,  namely, 
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1  5.  5  kcal/mole  for  reaction  (36),  and  32,0  kcal/mole  for  reaction  (37).  On 
the  other  hand,  the  dissociation  energy  of  Cl^  is  about  57  kcal/mole,  so 
that  the  energy  required  for  generation  of  each  mole  of  dissociated  Cl  atoms 
is  at  least  28.  5  keal.  Thus,  if  a  laser  is  used  as  the  energy  source  for  the 
generation  of  the  active  atoms  for  the  pumping  of  the  nonchain  chemical 
laser  reactions  of  the  type  (36)  and  (37),  the  required  laser  input  energy 
would  most  certainly  be  greater  than  the  resultant  c  hemical  laser  output 
energy.  Therefore,  one  can  quite,  safely  conclude  that  the  prospect  of 
achieving  pulse  energy  gain  from  any  single-phase,  nonchain  chemical 
laser  driven  by  another  laser  is  always  negative. 

d.  Prospect  of  Pulse  Energy  Gain  from  Chain  Reactions 

When  a  laser  is  used  as  an  initiation  energy  source  for  a  pulsed 
chemical  laser  based  on  chain  reactions  of  the  type 

F  +  H  -  HF  +  H  -t  31.7  kcal/mole  (38) 

H  +  l;z  -  I-IF  H  F  H  106.  7  kcal/mole  (39) 

or  their  isotopic  equivalent,  reactions  (1)  and  (2),  the  prospect  of  energy 
gain  appears  to  be  much  better  since  each  chain  center  atom  generated  by 
the  initiation  source  can  now  be  recycled  many  times  before  tlu:  chain 
terminates.  However,  in  pulsed  systems  where  the  reaction  produce  mole¬ 
cules  (e.  g.  ,  HF  or  DF  from  reactions  (38),  (39),  (1)  and  (2)  cannot  be 
physically  removed  from  the  lasing  region  within  the  pulse  period,  the 
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maximum  Laser  energy  output  that  can  be  extracted  through  the  consumption 
of  each  chain  center  atom  would  be  governed  by  the  collisional  de-excitation 
rate  of  tlie  upper  laser  state  as  well  as  by  the  reaction  chain  length.  We 
shall  examine  these  points  in  more  details  in  the  sub-sections  that  follow, 
since  these  points  are  quite  crucial  in  any  quantitative  prediction  of  single¬ 
phase  pulsed  chemical  laser  performance,  using  whatever  method  of 
initiation. 

(1)  Prediction  of  Pulsed  Chemical  Laser  Energy  Output  Based  on 
Chain  Reactions.  As  mentioned  in  the  introductory  section,  the  most  pene¬ 
trating  theoretical  analysis  of  single-phase  pulsed  chemical  laser  operation 
involving  chain  chemical  reactions  to  date  appears  to  be  that  of  Qraevskii,  ^  " 
even  though  the  analysis  was  limited  to  linear  systems  involving  very  few 
rate- limiting  reactions.  Oraevskii's  analysis  can  be  summarized  briefly 
as  follows: 

(a)  The  Single  Chain  Reaction 

An  idealized  three-level,  system  is  considered  here  (see  Fig.  1). 

Let  us  consider  the  following  chain  process, 
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where  and  B.,  are  the  initial  reactants,  AB(2)  is  a  molecule  in  the  upper 

laser  level,  AB(1)  is  a  molecule  in  the  lower  laser  level,  AB(O)  is  a  mole¬ 
cule  in  the  ground  state,  Q  is  the  energy  supplied  to  dissociate  one  oi'  the 
reactants  to  form  the  chain  centers,  and  M  is  any  one  of  the  reactants. 

And  k,  kg,  k^,  and  are  respectively  the  rate  constants  for  the  reactions 
(i),  (ii),  (iii)  and  (iv).  If  the  initiation  of  the  reaction  is  by  light,  then  k  is 
proportional  to  the  photon  density  and  the  pholod issociation  cross  section. 

The  rate  equations  describing  the  forward  reactions  of  these  processes  can 
be  written  down  for  the  concentrations  of  AB(2),  AB(1),  A,  B,  A ^  and 
B ^  ,  i.e.,  [AB(21],  [AB(1)],  [A],  [B],  (A^],  and  [B.J.  However,  [A and 
[B_J  are  assumed  constant  to  permit  a  linear  approximation  of  the  equations. 
Introducing  the  following  dimensionless  variables, 


|AB(2)] 

"  t 13  1 
2  0 


[  AB  (1 )] 

1  "  l»,L 


l»2] 

2  0 


0 1  , 


k  JB  i 


kjlM] 

kJliT 

2  2  Q 


k|A0j 

L.- 


the  rate  equations  can  be  written  as 


36 


dx . 


dx  . 


dT  V  ~  0 1X2 


dT  °]X2 


(4i: 


du 

dT 


-  v  -  au  +  Y 


dv 

dT 


=  -  v  i  a  u 


With  initial  conditions  such  as  x  =Xj=u=v  =  0  at  t=0j  and  the 
assumptions  that  the  excitation  short  pulse  duration  T  is  much  less  than 
the  relaxation  time  and  that  1  +  O  >:>  0  ,  the  following  expression  for  the 

maximum  value  of  population  inversion  is  obtained 


(x  -  x  ) 

^  1  max 


f° 

J  Y(T')dT' 


(42) 


where 


(1  -  In  2) 


a 


1  H  O 


is  the  chain  length.  It  is  obvious 


from  Eq.(42)  that  the  maximum  inversion  density  is  propo  "tional  to  the 

absorbed  number  of  pump  quanta.  Therefore,  even  for  long  effective  chain 

length  v  ,  =  t — ~ —  — .  the  output  power  of  the  chemical  laser  based  on  a 
eff  1  +  o  0  r  1 

simple  chain  reaction  will  be  determined  by  the  pump  intensity. 


(b)  The  Branched  Reaction 

From  chemical  kinetics  point  of  view,  the  branched  reactions  give 

grounds  for  supposing  that  the  number  of  active  particles  obtained  in  the 

final  stage  is  independent  of  the  energy  initiating  the  reaction  since  it  has 

37  , 

been  shown  for  such  reactions  that  the  rate  grows  exponentially  with  time. 
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To  obtain  a  population  inversion  in  the  system,  a  number  of  limitations  on 
the  rates  of  the  elementary  reactions  of  a  branched  chemical  process  are 
derived. 

First,  the  following  branching  scheme  is  assumed, 

k2 

(v)  AB(2)+RA  -  AB(1)+R+A 


where  RA  is  some  additional  reagent  that  dissociates  upon  collision  with 
AB(2),  By  incorporating  reaction  (v)  into  Eq.  (41),  the  following  system  of 
equations  are  obtained 


du 

dT 


=  v 


V-  (Oj  *».|  ^  , 

-  O'  u  I-  0  x  _  +  Y  , 


d  v 

dT 


-  V  i  OU 


(43) 


k'(RA] 

where  a  =  ~r — — ,  and  k  '  is  the  rate  constant  for  reaction  (v),  The 
2  k  1 13  |  2 

2  2  o 

solution  to  may  be  written  as 


x2  =  Ao 


ST 


(44) 


where  A  is  a  constant  which  depends  on  the  intensity  of  illumination,  and 
s  is  the  positive  root  of  the  equation 


3 

s  -f 


s 


2 


(1 


IfllOjlO  )  f  S  ( 1  IflK'Jjifl^) 


o  -  0 


(45) 
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and  it  follows  from  Eq.  (43)  that 


°1  +  a2  „  ST 

x,  = -  Ae 

1  s 


The  condition  of  exponential  growth  of  the  inverted  population,  i.e. 
(x^-  Xj)  >  0  implies  the  following  inequality 


s  >  °i  +  a2 


Therefore  a  general  condition  which  imposes  a  limitation  on  the  partial 
densities  of  the  reactants  and  the  rate  constants  of  the  reacting  mixture  is 
obtained. 

To  be  more  specific,  consider  the  following  case  of  low  branching 
rates  such  that  <<  o  ,  and  also  o  «  1.  From  Eq.(45),  we  obtain 


s  s  o  /  o  ^ 


Thus,  the  inequality  (47)  becomes 


aa_  >  a, 

L  X 


Equation  (49)  can  then  be  written  as 


(1  +  r,  4  ?.)2  k0k2 


where 


V) 


[RA]  ,  1  2J  • 

TbJ  ■  and  T|  =  T5J  ' 


Now  Eq.  (50)  implies  that  §  must  satisfy  the  following  inequality, 


where  ?  and  ?  are  the  roots  of  the  quadratic  equation 


? 2  +  2  ( 1  +  n  -  n  -2-1. 

\  2k, 


+  (H  h  )  -  0 


Since  physically  ?  >  0,  therefore  froin  (50) 


k  k'  ,  , . 

0  2  _  1  f  m 


Similar  conclusions  can  be  drawn  on  h  since  it  is  obvious  in  Eq,  (10)  that 
both  T|  and  ”  are  symmetrical. 

Now,  if  the  branching  scheme  is  the  following 


(vi)  A 13 (2)  +13  -  213  +AB(0) 


then  following  similar  analysis  as  before,  we  obtain  the  following: 


h  <  h  <  n 


where  P  and  p  are  the  roots  of  the  equation 
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n2  +  [a-^-l  „♦>=* 

l  h2  J 


and  we  must  have 


k0k2 


And  if  the  branching  of  the  reaction  is  due  to  the  dissociation  of  A, 


molecules : 


(vii)  AB{2)  -t  Az  =  2.A  +  AB(O) 


then  it  is  necessary  that 


koki 


T)  <  n  <  r| 


where  +1  and  h  are  the  roots  of  the  equation 


2  /,  \  , 
11  i1'— )+z 


r)  -i  1  +1=0 


It  can  be  observed  that  conditions  (5  3),  (55),  and  (57)  all  indicate 


that  in  order  to  have  exponential  growth  of  the  inverted  population  with  time 


41 


in  the  initial  stage  of  a  branched  chain  reaction,  the  effective  reaction  rate 

/  1  /2 

of  production  of  the  active  chain  centers  (k  k.,)  must:  be  greater  than 
the  relaxation  reaction  rate  (for  reaction  iv). 

A  numerical  solution  of  a  model  problem  is  solved  with  arbitrarily 
selected  rate  constants  demonstrating  the  independence  of  the  maximum 
inversion  on  the  intensity  of  initiation  of  reaction.  This  independence 
could  be  expected  from  the  exponential  growth  of  inversion  at  the  initial 
moments  of  time. 

(2)  Kinetics  of  the  HF  Chemical  Laser  System.  Experimental 

observations  of  laser  action  of  hydrogen  fluoride  v-is  first  reported 

3  8 

in  1967  by  Kornpa  and  Pimental.  They  observed  this  laser  action 

in  the  flash  photolysis  of  uranium  hexfluoridc/hydrogen  (deuterium) 

mixtures.  Since  then  several  other  fluoride  /hydrogen  mixtures  were  used 

39-43 

in  the  production  of  hydrogen  fluoride  laser  by  flash  photolysis.  The 

first  continuous-wave  HF  chemical  laser  was  reported  by  Spencer,  Jacobs, 


Mivels  and  Gross. 


They  used  a  supersonic  nitrogen  jet,  containing 


a  dilute  concentration  of  F  atoms,  flowing  into  an  ambient  hydrogen  atmos¬ 
phere  and  past  an  optical  cavity.  Molecular  hydrogen  was  then  injected 
through  perforated  tubes  transverse  to  the  jet.  Diffusion  of  H,  into  the 
jet  produces  HF  which  lases  in  the  optical  cavity. 


The  spectral  lines  observed  from  continuous  HF  and  DF  laser 
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outputs  are  given  in  Table  I  reproduced  from  Ref,  9.  From  the  intensities 
of  the  lines  shown  in  Ref.  9,  it  is  found  that  P^(5)  transition  at  2,  795  u 
arid  P^(6)  transition  at  2.  707  U  are  the  strongest,  in  7IF  laser.  It  is  also 
seen  that  all  the  lines  are  in  the  P  branches  of  the  various  bands.  This 
observation  is  confirmed  by  their  calculations  of  <he  sero-power  gam 
coefficients  of  the  transitions  in  a  given  band.  The  calculations  show  that 
the  P- branch  lines  generally  have  greater  gain  coefficients  and  the 
R-branch  coefficients  do  not  even  have  positive  values  unless  the  total 
inversion  in  the  populations  of  two  vibrational  states  is  quite  significant. 

The  proposed  chemical  chain  reaction  used  for  the  II,,  -  F  system 
proceeds  as  follows: 

(viii)  F  -f  M  2F  I  M  (Chain-initiation  step) 

(ix)  F  -1  H2  -  HI'  (v)  +  H  J 

|  (Chain-propagation  step) 

(x)  H  +  F  -  H.F  (v )  I  F  ) 

(xi)  HF(v)  H  M  ~  HF (v  -  1 )  I  M 

(xii)  HF  (v)  +  HF  (v ' )  =*  HF!v  +  l)HIF(v'-li 

(xiii)  II  H  H  +  U  *  +  M  \ 

O  I 

/  (Chain- te >  mination  step) 

(xiv)  F+F+M^F  +M  ) 

ti 

Reaction  (viii)  represents  the  dissociation  of  fluorine  molecule  by  collision 
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Table  I 


OBSERVED  SPECTRA  EROM  THE  CONTINUOUS  CHEMICAL  LASERS 


IIP 

DP 

Identification 

vibrational 

band 

Line 

Wave¬ 

length 

(u ) 

Identification 

vibrational 

band 

Line 

V  a  vo¬ 
le  ngth 
(u  ) 

1  -  0 

P(4) 

2.  640 

1-0 

P(8) 

3.  680 

P  (5 ) 

2.  67 3 

P  ( 9) 

3.  716 

P(6 ) 

2.  707 

P(1  0) 

3.  752 

P  ( 7 ) 

2.  744 

P(i  1 ) 

3.  790 

2  -  1 

P  (4) 

2.  76  0 

P  (1  2) 

3.  830 

P(5) 

2.  795 

2  -  1 

P(8 ) 

3.  800 

P(6) 

2.  832 

P(9) 

3,  838 

P(7) 

2.  871 

P(10) 

3.  876 

P(ll) 

3.  916 

P(1  2) 

3.  95  7 

3-2 

P  (8 ) 

3.  927 

P(9) 

3.  965 

P(10) 

4.  005 

P(1 1 ) 

4.  046 
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with  M,  a  collision  partner.  In  reactions  (ix)  to  (xii),  a  particular  vibra¬ 
tional  level  of  HF  is  indicated  by  v.  Reactions  (ix)  and  (x)  arc  the 
excitation  processes  of  HF  molecules;  reactions  (xi)  are  the  vibration- 
translation  (V-T)  transfer  equations  with  single  quantum  steps  only; 
vibration -vibration  (V-V)  exchange  of  a  single  quantum  is  represented  by 
reactions  (xii). 

Different  excitation  and  de-excitation  rates  for  the  H  -F  system 
44-47 

were  also  reported.  In  fact,  a  detailed  list  of  reactions  considered 

for  the  rate  equations  of  the  vibrational  levels  of  HF  molecule  together 

47 

with  selected  rate  constants  prepared  by  the  Aerospace  group  is  pre¬ 
sented  in  Appendix  1. 

(3)  "Inversion  Boundary"  for  the  Pulsed  H  -F  System.  Due  to 

L-.  Lj 

the  complexity  of  the  chemical  reactions  involved  in  the  HF  chemical  laser 
(H^-F^  system),  it  is  very  difficult  to  solve  the  resulting  non-linear 
chemical  rate  equations  without  resorting  to  numerical  techniques.  The 
simplified  treatment  of  Ref.  10  on  chemical  lasers  cannot  be  applied 
directly  to  the  HF  laser  without  making  much  more  elaborate  analyis  be¬ 
cause  the  HF  laser  is  not  a  simple  3-level  system  as  depicted  in  Fig.  1, 
and  many  more  reactions  are  involved  (Appendix  1).  Therefore,  the 
"inversion  boundary"  of  HF  laser  is  defined  here  as  the  dividing  line  in  the 
and  n  plane  which  indicates  the  possible  extent  as  to  how  the  IIF 
laser  may  progress.  (Here  n  denotes  the  number  of  type-C  molecule  or 
atom  per  unit  volume  at  any  given  instant.  ) 
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It  has  been  reported^  ^  ^  that  the  main  laser  actions  of  the 

HP'  chemical  laser  are  between  the  second  and  the  first  (and  also  between 

the  first  and  the  ground)  vibrational  states  of  the  HP  molecules.  For  the 

sake  of  simplicity,  let  us  only  consider  laser  action  between  the  second  and 

the  first  vibrational  states  of  the  HF  molecules.  Since  all  the  reaction  rates 

are  supposed  to  be  known  (Appendix  1),  the  detailed  rate  equation  for 

the  first  and  the  second  vibrational  state  population  can  be  written  down 

respectively.  These  two  equations  will  generally  involve  the  following 

unknowns:  the  various  vibrational  state  number  density  of  the  HP"  molecule, 

nTT„.  %  >  v  =  0,  1 ,  2 ,  3,  .  .  .  ,  the  hydrogen  atom  and  molecule  number  density 
PlF(v) 

n I j. ,  n  ,  respectively,  the  fluorine  atom  and  molecule  number  density 

.  ,  .  dnHF(2)  dnHF(v)  ,  , 

n  ,  n  respectively.  By  setting  - — - r-  - — — ,  an  algebraic 

J '  r  ^  "  cl  t  cl  t 

equation  results.  With  some  simplification,  this  algebraic  equation  can  be 

reduced  to  one  involving  only  nTTT,  and  n^_  .  This  equation  can  then  be 

H  P  r 

represented  by  a  curve  in  the  n__  -  n _ plane  using  n_  as  the  independent 

variable.  Such  a  curve  is  what  we  call  the  "inversion  boundary"  of  the 
system. 

(a)  Total  Inversion.  In  the  initial  stage  of  the  overall  reaction, 
reaction  (ix)  is  responsible  for  the  vibrational  excitation  of  the  HI"  mole¬ 
cule,  It  can  be  accordingly  expressed  as 


P  -I  11^  -•  HP’(v)  •!-  II  ,  where  v  ---  Q,  1,2,  3 


(GO) 


A  7 


The  laser  action  occurs  approximately  at  3  jj  (see  Table  I),  therefore,  it 
requires  about  9.  6  kcal/mole  to  excite  the  HF  molecules  from  their  existing 
vibrational  levels  to  the  next  higher  vibrational  levels.  The  heat  of  reaction 
for  (38)  or  (60)  can  excite  vibrational  levels  up  to  v  =  3  for  the  HF  mole¬ 
cule,  It  is  then  assumed  that  only  v  =  0,  1,2,  3  vibrational  levels  of  HF 
molecule  are  excited  in  the  following  formulation  of  the  rate  processes  for 
the  first  and  second  vibrational  levels  of  the  HF  molecule. 

How  the  rate  equations  for  the  first  and  second  vibrational  levels  of 
the  HF  molecule  can  be  w  Htten  down  according  to  the  reactions  listed  in 
Appendix  1. 


^  nHF  ( 1 ) 

^  =  k4(l)nFnH„  +k5(l)nHnF  +k7a(l)nHF(Z)nHF(0)  4k7a(2)nHF(2)nHF(>) 

Z  Z 


+  k 7b  ( 2 )  nHF  ( 2 )  nHF  ( 3 )  '  K7b  ( 1 )  nHF  ( 1 )  nHF  ( 2 )  '  k 7a  ( 1 )  ni IF  ( 1 )  nHF  ( 1 ) 


"  k7a  (2)  nHF  (1 )  niIF  (3)  "  k?c  (1 )  nHF  (1 )  nHF  (3 )  ~k6a(l )  ^niIF  (1 )  nHF 

4  nHF(l)nH'J  H  k6a(2)nHF(2)  (nHF  +nH  ’ 

z  ^ 

"  k6b(l)rHF(l)nF  +  k6b(2)nHF(2)nF  ”  k6c  (T  )  ^HF  (1 )  “f  2 
+  k6c(2)nHF(2)nF7  '  k6f(l)niIF(l)nH  +  k6f  (2)  nHF(2)”lI  (61) 

La 
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HF(2)  - 1  , ,  ,  , 

dt  ~k4(2)nFnH0  '  k5(2)nHnF  "  k7a  (1 )  UHF  (2)  nHF  (0)"  k7a  (2)  nHF  (2)  nHF(2) 

^  c. 

"  k7b  ( 2 )  nl IF  ( 2 )nHF  ( 3 )“  k 7b  (1 )  nHF  ( 1 )  nHF  (2 )  +  k 7a  ( I )  nHF  ( 1 )  nHF  ( 1 ) 


+  k7a(2)nHF(l)nHF(3)  +  k6a(2)  ^nHF(2)nHF  +  nHF(2)nH2^ 


ek6c(2)nHF(2)nF2  +  k6c(3)nHF(3)nF2  '  k6f  (2)  ni-IF  (2)  nH 


+  k6f(3)nHF(3)nH 


where  all  n, .  .  can  be  related  to  n 


HF(v)  HF  2L,  HF  (v) 


=  I 


by  the  following 


relationship, 


HF  (v) 


k4(\ 


By  using  the  values  of  rate  constants  in  lief.  7,  we  have 


nHF(0) 


U,  '1552  n. 


1IF(1) 


^  0.11 1  n,. 


HF  (2! 


0.555  n. 


HF  (3 ) 


°-2775  nHF 


The  atom  conservation  equations  for  the  hydrogen  and  fluorine  atoms 


can  be  written  as, 
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(65) 


”h  “  2nHi<t  =  0)  -  2nH2  '  "HF 


”F  =  2nF2(t  =  0|  '  2nF2"  "HF 


where  n  and  n  _  are  the  initial  number  densities  of  hydrogen 

hi  7(t  -  UJ  h  _  (t  -  U ) 

.  nH  ■  nF 

and  fluorine  molecules.  Since  - — - «  1,  and - 1,  we 

nH0(t  -  0)  n*,(t  =  0) 


have  approximately, 


nHF  +2nH2  “  2  nH2(t  =  0) 


"HF  2  nF  2  ^  °F  2  ft  =  0 ) 


Tlms’  U  I,H2(t  =  0)  "  nF2(t=  0)  '  thcn  nF2(t)  *  nH2(t)  •  Md  «»s..qucntly, 
nH(t)  "  nF(t)‘ 

With  the  initial  condition  that  nTT  .  ....  =  nT^  .  A,  ,  we  can  introduce 

H  (t  =0)  I‘  (t  =  0) 

the  following  dimensionless  quantities  into  Eqs.  (61),  (62),  (64),  (65),  (66): 


nF ^  (t  0)  nH2(t-0) 


nF  2  (t  =  0 )  nH2(t-0! 


nHF(v) 


v  n 


f2(o) 
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n 


H„ 


n. 


z  = 


nH  (t  =  0)  nF  (t  =  0) 
2  2 


Then  Eqs.  (61),  (62),  (65), and  (66)  become, 


dyl 

“dT=  k4(l)  XZ  +k5(D  XZ  +k7a(l)Yly0  +k7a(2)  Y2Y2  4  k7b(2)  Y2Y3 


k7b(l)  yly2  "  k7a(l)  Ylyl  “  k7a(2)  Y1Y3  '  k7c(l)  Y1Y3 


k6a(l)  (ylY  +ylz)  +k6a(2)  (y2Y  +  y2z)  “  k6b(l )  yl  X  +k6b(2)y2X 


k6c(D  ylZ  +  k6c(2)  V  “  k6f(l)  Y1X  +  K6f(2)  Y2X 


(69) 


dy. 


dt  ~  k4(2)  XZ  +k5(2)  XZ  "  k7a(l)  Y1Y0  k7a (2)  Y2Y 2  “  k7b (2)  Y 2Y : 


7b (2)  '  2  3 


"  k7b(l)  yly2  +k7a(l)  ylYl  +  k7a(2)  yly3  “  k6a(2)  Y2(y  +  z) 


+  k6a(3)  y 3 (y  +  Z)  ’  k6b(2)  Y2X  +  k6b(3)  y3X  '  k6c(2)  Y2Z 


+  k6c(3)y3a-k6f(2)yZX  +k6£(3)y3x 


(70) 


z  =  1  -  (x  +  y) 


(71) 


By  setting 


dy2  dyl 


dt  dt 


we  have 


k  xz  -  k  y“  -  k  f(y,  z)  =  0 
4/0 


(72) 
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v>  he  re 


k4  k4(2)  +  k5(2)  ~  k4(l)  ~  k4(2) 


k7  =  0.  59?  k  +  0.  5564  k 


7a(l)  T  k7a(2)  ~  0.  0308  +  Q.  308  k 


7b  (21 


kf,f  (y’ z)  =  y  [i.  11  k 


'6a  (2) 


-  0.  2775  k 


6a(3)  '  111  k6c(l)J 


+  yZlI'I‘(k6a(Z)  +  k6c(2,1-  0-  *775k6.(3,Jk«c(J,>-0.»»ft6jl(1)+k6cn)«] 
+yxfI-  U  (1'6b(2)+k6f(2)»-  111  <k6ia4k6m),-°-  2775  a'61a(3)+k6f(3>,J 


By  substituting  Eq.  (71)  into  Eq.  (72),  there  re: 


y2  <  *  2  k6bf  '  *W  *  1  2k( 


+  _  2  „  k4 

[2(k,7  +k/  )-k,  ]  x  "  2  -0 

‘  oa  bac  1 . 7  "'6a '  ~  ^bac J 


who  r< 


^  =  l-nk6a(2>-  °- 7775  k 


.  -  0.  Ill  k,  , 

6a(3)  6a.fi ) 


hac  "■  11  "W)-1  W-  °-  2775  lk6a(3.+k6c<3|l-  "•  H  Ik6aa)+k6c|l)] 


kM,t  “  11  lk/, 


6b(2)  6f(2)  ^  lkbb(3)+k6f(3)]  •  111  lk6h(1)  +  *6f(1)1 
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w  here 


k4  k4(2)  +  k5(2)  "  k4 (1 )  '  k4(2) 


k_  =  0.  597  k  +  0.  5564  k  -  0.  0308  k_  +  0.  308  k 

7  7a(l)  7a(2)  7c  ( 1 )  7b(2' 


k6£(y’7,)  =  y  t1'  11  '  °-  2775  kA,m  '  °-  111  Vwi  J 


6a  (2) 


6a  (3) 


6c  (1 ) 


+  yz[1*U  (k6a(2)+k6c(2))-°'  2775k6a(3)  +  k6c(3))-°-111  \a(l  )  +  k6c  (1  )V ] 
+yx  [1.  11  (k6b(2)  +  k6f(2))  '  °*  111  (k6b(l)  +  k6f(l))  '  °*  2775  (k6b(3)  +  k6f(3))^ 


By  substituting  Eq.  (71)  into  Eq.  (72),  there  results, 


(k  +  2  k  -  k  )  x  +  2  k 
2  ,  4  6bi  6ac  6ac 

y  4-  -  - - — - - y 

[2(k  +k,  )  -  k  ] 

L  7  oa  6ac  J 


k 


l2S+k6a»-k6ac! 


x2-2 


[aVVW1 


x  =  0 


whe  re 


(73) 


k6a  "  l-Uk6a(2,-  0-  2775  k6a(3,-0-nlk6a(l) 


k6ac  U  "WW  °'  2775  (3)  +  k6c  (3)k  11  ik6a(1  ,«6c!1 ,1 


k6bf  1,11  kf,b(2)  +  kfc£(2)k  °‘  2775  ‘k6b(3)+k6£(3)'  "  °‘  111  lk6b(l ) H  k6f(l)' 
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In  Eq.  (7  3),  all  the  rate  constants  k  are  known  functions  of  tempera¬ 


ture;  therefore,  once  the  temperature  of  the  system  is  given,  Eq.  (73)  can 

be  represented  by  a  curve  in  the  x-y  plane.  This  curve  is  the  "inversion 

boundary"  for  the  HF  chemical  laser  we  sought. 

In  Fig.  2,  the  HF  inversion  boundary  for  three  different  values  of 

the  translational  temperature,  T  =  300°K,  5  00°  K,  1200°K  are  presented. 

It  can  be  seen  from  the  figure  that  in  the  region  under  the  inversion 

boundary,  dnTIT,,_./'dt  is  greater  than  dn  ,  ,/dt  because  the  pumping 
H i*  (Z )  rix  ( 1 ) 

rate  by  n_  is  faster  than  the  de-excitation  rate  by  n  .  In  this  region, 

F  Hr 

population  inversion  in  the  system  is  most  likely  (i.  e.  ,  assuming  reason¬ 
ably  favorable  past  history).  Whereas  in  the  region  above  the  inversion 
boundary,  the  opposite  is  true.  As  expected,  the  inversion  boundary  moves 
upward  as  the  translational  temperature  increases  from  300°  K  to  SOO^K; 
but  the  upward  move  decreases  as  the  temperature  increases  further  from 
500°  K  to  1200°  K.  These  phenomena  are  because  the  excitation 
rate  constants  are.  much  steeper  functions  of  temperature  than  the  de¬ 
excitation  rate  constants  in  the  temperature  range  from  300  °K  to  500°  K. 

As  temperature  increases  further  from  500°  K  to  1200°  K,  the  temperatur e 
dependence  of  both  the  rate  constants  are  comparable,  hi  other  words,  as 
the  temperatur e  increases  from  500°K  to  1200°K,  the  de- excitation  rate 
constants  of  the  HF  chemical  laser  catch  up  slowly  with  the  excitation  rate 


constants,  hi  fact,  for  n„/nT_  .  A.  i  10  ,  it  is  seen  from  Fig.  2  that 

F  t7(t  =  0) 


Inversion  boundaries  as  functions  of  the  translational  temperature 
for  single-phase,  H^-Ft  chemical  laser  systems. 


the  inversion  boundary  at  1?00°K  is  actually  lower  than  the  inversion 
boundary  at  500° K. 

(b)  Partial  Livers  ion.  So  far,  we  have  only  discussed  the 
case  of  complete  inversion  (or  vibrational  inversion)  in  the  HF  chemical 
laser  system,  i.  e.  ,  n(v/)  and  n(v)  are  the  number  densities  of  the  upper 
and  lower  vibrational  states  v'  and  v  of  a  two-level  laser  system. 

In  fact,  it  is  not  always  necessary  to  have  complete  inversion  to 
ensure  laser  action  between  two  levels  of  the  same  vibrational  mode  in  a 

molecule.  Laser  action  may  be  obtained  also  if  there  is  only  partial 

48-50  ,  ,  ,  .  .  .  . 

inversion  (vibration- rotational  inversion),  1.  e.  , 

n(v'j')  -  n(v,j)  >  0  (7 


n(v ' )  -  n(v)  ^  0 


where  n(v,  j)  is  the  number  density  of  the  specified  vibrational  state  v 
and  rotational  state  j,  and  g  (j )  is  the  rotational  degeneracy,  2j  +1.  From 
statistical  mechanics,  we  know  that  if  the  rotational  states  are  in 


equilibrium  among  themselves  at  temperature  ,  n(v,  j)  would  be  given  by 


-1  -E0«)/kTr 

n(v,j)  =  n(v)  g(j)  Q^ot  e 
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The  rotational  partition  function  Q  and  the  energy  of  the  rotational 
state  j  ,  E^fj)  can  be  expressed  as 

kTr 

^rot  he  B 


E 0 (j )  -  he  B  j  (j  t  1 )  (relative  to  j  =  0  state  of  that  same 

vibrational  level  v) 


where  n(v)  =  2_,  n(v,  j)  is  the  total  numln  i  density  of  the  speci.iod  vibra- 

j 

tional  state  v,  lc  is  the  Boltzmann  constant,  h  is  the  Planck  constant,  c  is 
the  speed  of  light,  and  B  is  the  rotational  constant.  The  rotational  tem¬ 
perature  T  may  be  assumed  to  be  equal  to  the  translational  temperature  T. 

Now  in  order  to  find  the  "partial -inversion  boundary,  "  we  have  to 
use  the  following  condition: 

dn(v  ' g(j  ')  dn(v,  j) 

dt  'g(j)  dt  U  '} 


Assuming  that  the  temperature  is  independent,  of  time,  we  differentiate 
Eq.  (75)  with  respect  to  time  for  both  n(v/,j/)  and  n(v,j) 


dn(v/,  j') 
dt 


dn  (v ' ) 
dt 


^  ^o't1 


-E0(j/)/kTr 


(77) 


dn  (v,  j )  dn(v )  ,M  ^  - 1 

~dl~“  =  “d“^(j)  °rot  ° 


(78) 


By  substituting  Eqs.  (77) 


and  (78)  into  Eq.  (76)  }  we  obtain  the  following 
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An  algebraic  equation  similar  to  Eq.  (73)  can  be  obtained  with  different 
coefficients  by  substituting  Eqs.  (69)  and  (70)  into  Eq.  (BO). 

The  "partial- inversion"  boundary  for  t>e  HE  chemical  laser  are 
calculated  for  three  different  temperatures,  T  =  300°  K,  500°K  and  1200°  K, 
for  comparison  with  the  "complete-inversion"  boundary.  These  are  the 
dotted  curves  in  Fig.  2.  It  is  seen  that  the  "partial-inversion"  boundary  is 
generally  higher  than  the  "complete-inversion"  boundary.  That  is  to 
say,  under  identical  conditions,  it  is  easier  to  obtain  partial  population 


inversion  than  complete  population  inversion.  But  the  increase  is  not  very 

substantial;  at  T  -  300°K,  the  increase  is  about  a  factor  of  2,  at  T  =  500°K, 

the  increase  is  about  a  factor  of  1.3,  and  at  T  =  1  200°K,  the  increase  is 

only  about  a  factor  of  1.1.  Actually,  this  is  obvious  from  Eq.  (80)  because 

—  3  0  7  /  X 

of  the  Boltzmann  factor  e  .  As  the  temperature  increases,  the  effect  of 

the  Boltzmann  factor  drops  off  exponentially.  In  other  words,  at  high 
temperatures ,  complete  inversion  is  as  easy  to  obtain  as  partial  inversion. 

But  the  effect  is  much  more  pronounced  at  low  temperatures  where  partial 
inversion  is  much  easier  to  obtain  than  the  comolctc  inversion. 

(4)  Estimate  of  Maximum  Energy  Gain  Potential  for  Pulsed  H^-F^ 

I— cr  Systems.  J?hc  ''inversion  boundary  ol  the  Ill  chemical  laser  ]u s t 
presented  can  be  related  to  the  maximum  efficiency  or  energy  gain  potential  of 
the  laser  system  in  the  following  way.  The  inversion  density  An  is  defined  as 


An 


nIIF(2)  '  nlIF(l)  ’  nHF(2)  <  n ITF 


(81) 


where  ^  n,,,,  ,  is  the  total  number  density  of  the  T II*  molecule. 

HI*  '  Hr  (v) 

v 


Now  the  total  pulse  energy  up  to  time  t  it 


e  (t )  =  hv  j  dt 

„  dt 


-v  hv  dn 


HE 


(82) 


T  ^  as  the  time  when  the  gain  of  the  laser  system  is 
h? 


5  H 


We  want  to  define 


zero,  i.  e.  .  t  =  T  when  An  =  0,  and  t  as  the  time  when  — ■ =  0. 

g  =  0  1  dt 

Therefore,  the  total  pulse  energy  of  the  laser  is 


E!t  )  =hv  I 
B  0 


g  0  d  (An)  fg° 

~dT  dt  <  hV  J  dnHF  =  hV  nHF(Vo: 


(83) 


nHF^'l^  can  re*ated  to  nHF(T  _  in  the  following  way, 

g 


n 


HF  (T . 


n 


HF(t  ) 
g  =  0 


(84) 


Since  is  the  tune  when  d  (An ) /dt  =  0,  it  is  reasonable  to  assume  that 


T  is  proportional  to  the  excitation  time  of  HF(2),  T 
1  cx 


And 


T  =  V  +  T  where  T  is  the  time  duration  starting  from  d(An)/dt  =  0 

^  —  o  1  4  C* 

until  An  =  0;  thus  *s  assumed  to  be  proportional  to  the  fastest  relaxa¬ 

tion  time  q  HF(2),  T  , 

rel 


n 


HF  (2) 


1  ex  k4(?.)  nrnn 


(85) 


.  _  t  ^ _ 1 

2  rel  k.  , n 


6a  (2)  II, 


Thus. 


59 


HF  (t 


HF(T1 


-2-  +  1  +  - - 

T  ,  T 

1  2 


1  ,  k4(2)nF 

k6a(2)  nHF(2j 


1+1.8 


k4(2)  nF 

hlUzi  ^ 


Now  we  want  to  deiine  the  maximum  laser  efficiency  or  energy  gain 


potential  as 


laser  output 


HFV-o’ 


max  max,  chem.  energy  input  3  nT,  ,,  . 

1  .  ( 0 ) 


i  nHF(Tg-o’  nRF,Ti 


3  “iIFfTji  nl'2(r; 


when;  nIT,  , .  v/a_  ....  as  a  function  of  ,n_  /  n„  ...  .  is  the  curve  in  the 
Hl-('Tj)  F  1<  v(0) 

"inversion  boundary  ’  plot  of  Fig.  2.  Therefore,  the  maximum 

laser  efficiency  ti  is  directly  proportional  to  nrT_.  ./nT,  ,,,  . 

7  max  ’  1  +IF (+  ^ )  I-  (0) 

For  most  of  the  operational  HF  lasers,  me  ratio  of  ihe  free  fluorine 
atom  concentration  to  the  hydrogen  molecule  concentration  is  about  4%, 


i,  e.  n  /n.  ....  4%.  From  the  rate  constants  compiled  in  Append*  2 

I  i  2 (0) 

and  the  "inversion  boundary"  plot,  it  can  be  shown  that  fer  n^., !x\y  4%, 
at  T  =  300  °K 
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nHF<Tg  =  o) 

DHF(ti) 


2.  8 


max 


nHF(r  ) 

0.9 - - 

nFz(0) 


at  T  =  5  00°K 


nHF(T  ) 

_ g  =  0 

nHF(T1) 


11 


max 


nHF(Ti) 

JV2(o) 


at  T  =  1 200  °K 


nHF(x  ) 
8  =  0 

nHF  (t  1 ) 


n 


HF(T1) 


max 


n 


f2(0) 


Thus,  the  ’’inversion  boundary"  plot  can  serve  as  an  approximate  maximum 

laser  efficiency  plot  as  a  function  of  n^/n^  ^  . 

3 

5.  LASER  IGNITION  OF  TWO- PHASE  SYSTEMS 

a.  Definition  of  Two-phase  Pulsed  Chemical  Laser  Systems 

For  the  purpose  of  the  present  study,  wo  define  a  two-phase  pulsed 
chemical  laser  system  as  a  pulsed  chemical  laser  system  in  which  the 
lasing  medium  is  made  up  of  a  reactive  gas  mixture  plus  some  finely- 
dispersed,  submicron  size  solid  particles  or  liquid  droplets.  The  apparent 
advantages  of  two-phase  systems  over  the  commonly  considered  single¬ 
phase  systems  (i,  e.  ,  gaseous  mixtures  only)  are  many: 
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(1)  The  absorption  of  infrared  and  visible  radiation  by  small  particles  is 
generally  broad-band,  which  allows  a  wide  selection  of  the  chemical 
and  surface  properties  of  the  absorber. 

(2)  By  suitable  choice  of  particle  size  and  surface  properties,  the  heated 
particles  may  be  used  to  provide  ?.  large  catalytic  surface  area  for 
chain  initiation  and  for  subsequent  chain  branching  without  the  need 

for  grossly  heating  the  entire  gas  mixture  during  the  initiation  process. 

(3)  By  suitable  choice  of  chemical  properties,  the  heat  of  reaction  from 
combustible  particles  may  be  utilized  for  generation  of  a  high  concen¬ 
tration  of  active  atoms  essential  for  efficient  pumping  of  most  known 
chemical  laser  systems. 

Point  (3)  is  especially  attractive  from  the  point  of  view  of  potential 
laser  energy  amplification  due  to  the  very  large  ratio  between  the  heat  of 
combustion  of  most  metals  and  the  energy  required  to  heat  the  same  metals 
to  their  cor  responding  ignition  temperature,  starting  from  room  temperature 
condition.  Point  (1)  is  helpful  when  infrared  lasers  are  to  be  used  as  initia¬ 
tion  sources  in  view  of  the  generally  narrow-band  nature  of  the  absorption 
process.  (For  example,  the  cascade  excitation  dissociation  process  reported 
in  Reference  34  was  a  resonant  absorption  process,  and  hence  required  acci¬ 
dental  matching  of  the  incident  laser  frequency  and  the  molecular  vibrational 
frequency.  )  Point  (2)  could  be  useful  if  long  cnains  could  be  utilized  for 
chemical  laser  pumping  and  if  the  rate  ox  chain  initiation  in  the  otherwise 


purely  gas -phase  reaction  wore  too  slow. 
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We  shall  demonstrate  these  points  in  the  subsections  that  follow, 
b.  Optical  Properties  of  Two-phase  Systems 

The  absorption  and  scattering  properties  of  finely-dispersed,  sub¬ 
micron  size  particles  (or  liquid  droplets)  at  infrared  frequencies  can  be 

calculated  according  to  classical  electromagnetic  theory.  According  to  such 
52 

theory,  the  total  extinction  cross-section  Q  t  for  a  particle  of  equivalent 

3 

spherical  volume  V  =  4rra  /3  and  with  mean  radius  a  smaller  than  the 
radian  wavelength  X  =  \/2TT-l/k  of  the  incident  radiation  can  be  repre¬ 
sented  by  an  asymptotic  series: 


where  Re  and  Im  denote  taking  real  and  imaginary  part,  respectively, 
x  =  ka  =  a/X  is  the  radius-to-wavelength  ratio,  and  m  is  the  complex  index 
of  refraction  of  the  material  that  makes  up  the  particle.  The  first  term  is 
the  main  absorption  term,  and  the  second  term  is  the  main  scattering  term 
which  gives  rise  to  the  Rayleigh  scattering  formula  when  m  is  a  real  constant. 
By  writing  m  =  n.  -  itt  ,  where  n.  denotes  the  ordinary  index  of  refraction 
and  K  denotes  the  bulk  extinction  coefficient  of  the  material,  one  obtains 
from  Eq.  (88)  the  following  expressions  for  the  absorption  and  scattering 


cross  sections  of  the  particle  in  the  long  wave  limit  (i.  e.  ,  ^  »  a,  or 
x  «  1). 


Q 


abs 


24  ?vk: 


iia 


2  2  ,2  22 

(A  -  K  4  2)  4  4  It 


ka 


(89) 


Q 


scat 


f  (n  *  ) 


7  7  7  >  7  7 

rra  [(a  -  K  +  2)“  4  4«  J  K  ] 


—  (ka) 


(90) 


where 


8  2  2  4  22  2  22 

f(n,K)  =—  (n'  +  *  )  +  Z{n.  -  K  )  [ ( +v  )  -2] 

-  3  -  46  K 2  -  3  ^  4  4 


If  the  particle  material  is  a  low-loss  dielectric  (o.  g.  ,  quarts,  glass,  ionic 
crystal,  water,  etc.  ),  the  extinction  coefficient  is  generally  much  smaller 
than  the  index  of  refraction  (i.  e.  ,  K  «  >n).  For  such  particles,  the  two 
cross  sections  become  approximately, 


Q 


abs 


tt  a 


24  n.  K 
(a"  4  2)‘ 


--  ka 


(91) 


q  8  9  4  2 

scat  8 ( n  H  Zn  -  3  -  4 n.  +4;  ,4 

=“  - - - -  (ka) 


tt  a 


2 


(92) 


1  (A"  +  2) 


For  metallic  particles,  A  and  K  are  nearly  identically  given  by 

1/2 


A-  ^ 


a  X 
c 


(93) 
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where  o  is  the  electrical  conductivity  of  the  metal  expressed  in  Guassian 
units,  X  is  the  wavelength,  anl  c  is  the  velocity  of  light  in  vacuum. 

(Note  that  one  M.  K.  S.  unit  of  conductivity,  i.  e.  ,  1  mho /meter,  is  equivalent 

-  9  -1 

to  9  X  j-0  Gaussian  units  of  conductivity,  which  has  the  dimension  of  sec  .  ) 

For  a  typical  good  conductor,  such  as  aluminum,  the  numerical  value  of  a 

17  - 1 

at  room  temperature  is  typically  3  X  10  sec  .  Thus,  at  a  typical  infrared 

frequency  f  =  c/X  -  3  X  1 0 * sec  ^  (corresponding  to  X  -  10  p),  the  numerical 

2 

value  of  rt-  and  K  for  a  good  conductor  is  of  the  order  of  10  ",  which  is  much 
greater  than  unity.  Accordingly,  the  absorption  and  scattering  cross  sections 
for  submicron  size  metallic  particles  at  infrared  frequencies  ate  approxi¬ 
mately  given  by 


ahs 

2 

TT  3, 


6  C  . 

— T  K3. 

U  K 


(94) 


Q 

s 


cat 


2 

n  a 


(95) 


It  is  interesting  to  note  that,  since  Q  is  proportional  to  the  first 

cii  D  G 

power  of  ka  while  Q  ,  is  proportional  to  the  fourth  power  of  ka,  the 
scattering  cross  section  decreases  much  more  rapidly  with  particle  size 
than  the  a.bsorption  cross  section,  so  that  at  sufficiently  small  values  of  ka, 
the  scattering  cross  section  will  become  insignificant  in  comparison  with  the 
absorption  cross  section.  In  other  words,  a  two  -phase  mixture  composed 
of  very  fine  suspended  particles  will  essentially  behave  like  a  good  absorbing 
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medium  rather  than  a  strong  scattering  medium  at  infrared  frequencies. 


From  Eqs.  (91)  and  (92),  it  is  seen  that  for  dielectric  particles,  the  scatter¬ 
ing  cross  section  begins  to  fall  below  the  absorption  cross  section  when 
1/3 

ka  <  3  '(  •  (Note  that  m  -1,4  for  most  dielectric  material  at  infrared 

frequencies.  )  Similarly,  from  Eqs.  (94)  and  (95),  one  may  deduce  that  the 
scattering  cross  section  for  metallic  particles  will  fall  below  the  absorption 
cross  section  when  ka  <  1.  3  (c/ct  X  )  . 

Another  interesting  point  to  note  from  Eqs.  (91)  and  (94)  is  that,  if 
we  divide  these  equations  through  by  the  mean  particle  radius  a,  one  obtains 


Q  . 
abs 


48  tt  ■n'r 

?  ? 
(n.  H  2)  X 


1  2  IT  c 


a  \ 


for  dielectric 


for  metal 


(96) 


which  states  that  the  absorption  cross  section  for  very  small  particles  is 
proportional  to  the  volume  of  the  particle  at  a  given  frequency,  no  matter 

•*  -  -T— >  — -  4-  1  —  2  u.  d  -f  1  'v  4-  ci  v  ^  1 

wlid-t-  j-ii.Fi  wi  title  j.  x  icx  1,0  j.  i.ttJ-1 


Using  handbook  values,  the  extinction  coefficient  K  and  index  of 
refraction  x-  are  plotted  in  Figs.  3  and  4  as  functions  of  the  incident 
electromagnetic  wavelength  X  (in  vacuum)  over  the  infrared  region  for 
some  ionic  crystals  and  metals.  From  these  values  of  and  k  ,  the 
absorption  cross  section  for  very  small  size  particles  made  up  of  these 
mate  rials  (i.e.,  ka  «  1)  has  been  calculated  according  to  Eq.  (96)  as  a 
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X  ,  MICRON 


Extinction  coefficient  fC  ds  a.  function  of  wavelength  for  some 
ionic  crystals.  The  index  of  refraction  n  for  these  crystals 
has  a  nearly  constant  value  of  1,  4  in  the  infrared  range  of 
frequency. 


function  of  incident  wavelength.  The  results  are  plotted  in  Fig.  5.  From 
these  values  of  the  absorption  cross  section,  one  can  readily  calculate  the 
e-folding  absorption  length  scale  for  a  monochromatic  beam  of  electro¬ 
magnetic  radiation  propagating  through  a  two -phase  mixture  consisting  of  a 
transparent  gas-phase  component,  and  a  finely-dispersed  solid-phase 
component  of  n  particles  per  unit  volume, 


Since  the  volume  fraction  of  solid  particles,  0  ,  in  such  a  two-phase  mix- 

P 

ture  is  simply  given  by  the  product  of  n^  and  the  volume  per  particle 

3 

V  ~  4TTa  /3,  and  the  absorption  cross  section  per  particle  is  proportional 
to  the  particle  volume  according  to  Eq.  (9b),  one  can  subsequently  calculate 


the  value  of  0  for  yielding  a  desired  absorption  length  scale,  such  that 

P 


?  2 
(rc  +  2)  X 

36  rr  A  K 


,  T  -  4TTa 

p  abs  3  Q 

abs 


for  dielectric 


for  metal 


The  value  of  0  so  determined  for  an  absorption  length  scale  L  ,  =10 

p  abs 

meter  in  two-phase  mixtures  containing  the  indicated  solid  particles  is 

plotted  in  Fig.  6.  Since  the  bulk  density  of  most  solids  is  typically 
- 1  3 

10  mole/cm  while  the  standard  density  for  any  perfect:  gas  at  room  tem- 

-  5  3 

perature  and  atmospheric  pressure  is  about  4  X  10  mole/cm  ,  the  curves 
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X,  MICRON 


Absorption  cross  section  as  a  function  of  wavelength  for  sonic 
submicron  size  particles  (ka  <<  1)  of  some  ionic  crystals, 
metals,  and  graphite. 


\  MICRON 

Particle  volume  fraction  in  a  two-phase  mixture  which  would  yield 
an  e-  folding  absorption  length  scale  of  10  met  ers  at  room  tempera 
ture  for  some  ionic  crys'als,  metals,  and  graphite  as  a  function  ul 
wa\ elength,  neglecting  gas-phase  absorption  and  assuming  ka  <<  1 
n  i/.  e  na  r  t  it  1  es . 


rs^rs^ssaemm 


.shown  in  l'ig .  6  indicate  that  at  the  CO.,  lasev  wavelength  of  10.6  U  ,  the 
mole  fraction  of  solid  particles  required  to  yield  an  absorption  length  scale 
ol  10  meter  in  a  t-vo- phase  mixture  at  standard  atmospheric  pressure  and 
temperature  is  about  unity  when  the  particles  are  made  ol  highly  conducting 
metals,  such  as  aluminum  and  magnesium. 

c.  Solid  Particles  as  biert  Catalys'  for  Chain  Initiation  and  Chain 
Branching 

It  is  well  known  chat  hot  solid  surface  a  can  act  as  catalytic  surfaces 
for  dissociation  of  gaseous  molecules.  Dissociation  of  II  by  heated  tungsten 
is  just  one  of  the  examples.  In  gcr.eral,  one  may  expect  that  the  rate  of  dis¬ 
sociation  or  catalytic  efficiency  would  be  a  sensitive  function  of  tlm  solid 
surface  temperature,  as  observed  :n  the  dissociation  of  '■  i  by  tungsten. 

Thus,  the  large  surface  area  provided  by  suddenly  heated,  finely -dispersed 
solid  particles  in  a  two-phase  mixture  can  conceivably  he  utilized  as  a 
catalyst,  for  <.  hemical  chain  initiation  and  branching  in  the  gas -phase  compo¬ 
nent.  Of  course,  in  order  to  utilize  such  catalytic  effect  in  a  pulsed  chemical 
laser  system,  the  solid  surface  must  be  initially  passive  and  inert  so  that  the 
two-phase  mixture  would  be  stable  before  initiation. 

Detailed  calculation  of  the  surface-catalized  chain  initiation  and  chain 
branching  rate  would  require  a  specification  of  the  particle  heating  rate  as 
well  as  adequate  knowledge  about  the  temperature-dependent  catalytic  prop¬ 
erties  of  the  surface  of  interest.  While  such  catalytic  properties  may  be 
diflbult  to  get,  especially  for  particular  combinations  of  gases  and  sui  faces 


that  have  not  yet  been  studied.  The  initial  rate  of  heating  of  the  solid  particles 

by  a  suddenly  turned  on  beam  of  electromagnetic  radiation,  on  the  other  hand, 

is  relatively  easy  to  calculate.  For  example,  by  neglecting  scattering, 

re-emission,  and  conductive  heat  loss,  one  can  readily  write  down  the 

energy  balance  equation  for  the  small  solid  particles  of  individual  volume 

V  =  4naJ/3  and  of  buln  specific  heat  per  unit  volume  PC  during  the 

v 

passage  of  a  laser  beam  of  local  intensity  I, 


0  C  v  ^  4  IQ 

v  dt  abs 


Thus,  f.l.e  rate  of  tompcraiurc  rise  per  unit  beam  intensity  is  simply  given 
by 


1  dT  _  Qabs  _ _ 3__  Qab.; 

I  dt’  "  pc  V  4  pc  ’  "3 

v  v  n  a 


(100) 


The  inirial  particle  heating  1  ate  so  calculated  tor  the  various  particle 

material  and  absorption  cross  section  shown  in  Fig.  3  is  plotted  here  in 

Fig .  7.  it  is  seen  that  for  10.6  w  radiation,  fcne  heating  rate  for  aluminum 

6  2 

particles  by  a  beam  of  10  watts /cm  would  be  about  FK/u^ec,  whereas  for 
graphite,  the  corresponding  heating  rate  would  be  300 °K /usee.  At  solid- 
state  laser  frequencies  (.'■  ~  1  jj  ),  the  heating  rate  would  be  about  50  times 
faster, 

A  voter  the  chain  center  atoms  are  generated  by  the  heated  solid 
particles  through  surface  catalysis,  the  problem  of  subsequent  chain  center 


7  3 


diffusion  into  the  gas-phase  mixture  which  fills  up  the  interparticle  space 

can  be  treated  as  a  problem  of  spherical  diffusion  involving  conserved 

passive  scalar  quantities  (since  the  chain  center  atoms  arc  simply  being 

recycled  during  the  chain  propagation  step  in  nonbranching,  non- 

54 

terminating  chain  reactions):  It  can  be  shown  that  the  three-dimensional 
spherically  symmetric  diffusion  equation  with  constant  diffusion  coefficient 
may  be  written  in  the  following  way. 


(101) 


where  r  is  the  radial  distance,  n  is  the  concentration,  and  D  is  the 
binary  diffusion  coefficient. 

Consider  now  some  small  spherical  particles  distributed  randomly 

in  a  medium  of  dctonable  gas  mixture  such  as  IF /F  .  When  the  surface 

of  any  one  of  the  particles  is  uniformly  heated  to  some  temperature 

instantaneously,  dissociation  of  F  will  take  plat  e  on  the  surface  at  some 

value  of  T  .  Equation  (101)  with  some  appropriate  initial  and  boundary 
a 

conditions  may  be  used  to  describe  approximately  the  diffusion  process  both 
in  time  and  space  of  the  fluorine  atom  concentration  on  the  particle  surface. 
This  approximation  is  based  on  the  "equivalent  sphere"  concept  which 
assumes  that  spherical  symmetry  is  valid  in  this  case  up  to  a  certain  point 
in  space  for  each  particle,  A  schematic  diagram  of  the  problem  is 
illustrated  in  Fig.  8. 

Assuming  that  the  concentration  of  the  fluorine  atom  on  any  spherical 
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a/xw-si  * 


surface  remains  constant,  and  that  there  is  no  flux  at  the  imaginary  outer 
spherical  boundary,  we  can  set  up  the  proper  boundary  conditions  as 
follows. 


nF(t,  r) 


d_ 

Or 


nF(t, 


n 


Fa 


r)  =  0 


at  r  =  a  for  all  t 


at  r  -  t  +a  for  all  t 


(HU) 


And  the  initial  condition  that  there,  is  no  fluorine  atom  concentration  any¬ 
where  in  space  at  t-0  can  be  represented  as 

n_(t,  n)  =  0  at  t  -•  0  for  all  r  (103) 

F 

Similar  solution  of  the  problem  can  be  obtained  if  the  following 
dimensionless  quantities  are  introduced, 


(5  -  t) 


(r .  t) 

V 


(10d) 


By  substituting  the  expression  in  Eq.  (104)  into  Eqs.  HrU),  (102)  and  (103),  we 
obtain  the  following  uondimeiisional  equation  and  Us  initial  and  boundary 
conditions. 


d_ 

Ijt 


4 


(105) 


v  ith 


1.  C.  (t,  5  )  -  0  at  t  =  0  for  all  5  (106) 

B.C.  (T  +  ~)  rip  =  —  at  =  ^  f°r  a1^  t 

(107) 

9n 

( ?  +  t)  -  0  at  ?,  =  1  for  all  t 

v  st!  a-p 

This  may  be  transformed  into  a  problem  involving  diffusion  in  a  plate  of 
unit  thickness  (the  "equivalent  plate")  under  specific  boundary  conditions  by 
making  the  substitution  u  =  (?  +  V  )  np  *  Then  kqs.  (10<),  (105),  and  (106) 
become 

) 

du  .  (10b  i 

at  a? 2 

\v  ith 

I.  C.  u(t,  5)  =  0  for  0  <  ?  <  1  at  t  -  0  \109) 

B.  C.  u(t,  0  )  =  a/£  at  5-0  lor  all  t 

(HO) 

~  -  — — —  -  0  at  5  =  l  lor  all  1. 

°r-  [*  +  *\ 

Notice  here  that  the  second  boundary  condition  that  there  is  no  flux  at  the 
outer  boundary  in  the  original  problem  becomes  a  radiation  boundary  condi¬ 
tion  (in  the  sense  of  Ref.  54)  after  the  transformation.  According  to  R.ef,  54, 
the  solution  to  the  problem  represented  by  Nqs.  (10b),  (109),  and  (lib)  is 


i  i/i  ^  V  2(32  +12)  sin(j3  5  )  -02t 

u  1  +  L(1  -  ^  }  \  n  n  n 

a/i  1  +  L  n  =  1  P  (L  +L  +  p  ) 

n  n 


dll) 


where 


L  =  - 


1 

I 


and  6  are  the  positive  roots  of 
n 


3  cot  3  +  L  =  0 


(112) 


Therefore,  the  fluorine  atom  concentration  is  obtained, 


nF  (: ,  t) 


,  V  -  2(3  2  +  L“)  sin(g  5)  -p  2t 

a/i  >  n  n  n 

a  t~4  2  2  e 

§  +  7  n  =  l  P  (L+I/-+3*) 

f  n  n 


(113) 


Numerical  solutions  of  n  '  for  f /a  =  10,  100  with  different  values 
of  t  are  shown  in  Figs.  9  and  10.  It  should  be  remarked  here  that  for 
large  values  of  t,  the  infinite  series  in  Eq.(H3)  converges  very  rapidly. 
But  for  small  values  of  t,  more  and  more  terms  in  the  series  have  to  be 
used;  in  our  case,  for  t  =  0.  01,  six  terms  of  the  series  are  used. 
Fortunately,  for  even  smaller  values  of  t,  we  can  approximate  this  diffu¬ 
sion  problem  with  finite  boundaries  by  the  semi-infinite  diffusion  problem. 
Physically,  for  the  case  of  small  values  of  t,  the  concentration  does  not 


have  time  to  diffuse  out  far  enough  to  "feel"  the  presence  of  the  outer 
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Normalized  solution  for  transient  spherical  diffusion  of  cha in 
center  atoms  from  a  catalytic  particle  surface  located  at 
radius  a  into  a  gaseous  volume  of  interparticle  separation 
distance  2  9  ,  where  !.  / a  ~  10.  Other  nomenclature: 

°  F 

n  .  =  - - —  =  normalized  chain  center  concentration 

i  (n  Tr) 

1  a 

r  —  cX 

f  =  -  s  normalized  radial  coordinate 

t 

=  normalized  time  (D  s  diffusion  coefficient) 


t 


Fig.  10  Normalized  solution  for  transient  spherical  diffusion  of  chain 
center  atoms  from  a  catalytic  particle  surface  located  at 
radius  a  into  a  gaseous  volume  of  interparticle  separation 
distance  2f  ,  where  £  /a  -  100.  Other  nomenclature: 

nF 

nw  =  - - r—  £  normalized  chain  center  concentration 

1'  (n  ) 
r  a 

3^  “  3. 

%  2  -  ^  normalized  radial  coordinate 

*  £ 

— - -  '=  normalized  time  (  11=  diffusion  coefficient), 

£  /  D 
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t 


boundary.  This  is  verified  by  taking  the  solution  of  the  semi- infinite 
diffusion  problem  at  t  =  0.  01  to  compare  with  the  solution  given  by 
Eq.(ll3),  And  the  comparison  shows  very  little  difference  between  the 
two  solutions.  Thus,  for  t  «  1,  we  have 


a 


where  erfc(x)  is  the  complem  tary  error  function  of  argument  x, 
d.  Two-phase  Systems  Involving  Solid  Particle  Combustion 

(1)  Energy  Gain  Potential.  When  the  heated  particles  in  a  two- 
phase  chemical  laser  system  is  used  simply  as  an  inert  catalyst,  the 
energy  gain  potential  of  the  two-phase  system  is  not  very  different  trom 
that  of  a  single-phase  system  using  the  same  gas-phase  reactive  mixture. 
The  reason  is  that  the  energy  required  to  generate  the  dissociated  atoms 
must  still  be  supplied  by  the  initiation  source,  so  that  the  pulsed  energy 
gain  potential  would  still  be  governed  by  collisional  dc-excitation  and 
useable  chain  length  considerations  as  discussed  in  subsections  3.  c  and  3.  d. 

If  the  solid  particles  suspended  in  the  two-phase  system  can  be  made 
to  burn  and  the  subsequent  release  of  combustion  energy  could  be  utilized 
for  further  chemical  pumping  of  the  laser  reactions,  the  prospect,  of 
energy  gain  appears  much  better.  The  improved  p rospect  of  energy  gain 
derives  from  the  fact  that  the  heat  of  combustion  of  most  metals  and  some 
solid  compounds  is  often  much  greater  than  the  energy  required  to  heat  the 


same  mas  of  metal  or  compound  to  their  respective  ignition  temperature, 

starting  from  room  temperature.  This  is  illustrated  in  Table  2  in  which 

some  known  ignition  temperature  T^  for  various  solid  dust  particles  in 

air,  the  heating  energy  per  mole  AE^  for  raising  the  solid  temperature 

from  T.,  =  300°K  to  T  ,  and  the  heat  of  combustion  II  per  mole  of  the 
u  c  c 

■  •  it-  56 

initial  solid  in  O  I\,  and  Cl  are  listed.  Even  though  the  ignition 

£>  t-.  Ci 

temperature  for  these  solids  in  F-,  and  Cl.,  are  not  known,  they  are  likely  to 
be  lower  than  those  corresponding  to  the  air  case.  Thus,  the  heat  of  com¬ 
bustion  of  the  more  energetic  solid  fuel  such  as  aluminum  and  boron  in 
fluorine  is  generally  greater  than  the  energy  required  to  bring  the  solid 
particles  up  to  ignition  temperature  by  a  factor  of  about  100.  For  combustion 
in  chlorine,  the  energy  ratio  is  of  the  order  of  30, 

The  dissociation  energy  for  F and  CI^,  is  about  31  and  57  kcal/mole, 
respectively.  Thus,  the  heat  of  combustion  from  1  mole  (or  27  grams)  of  Al 
in  an  F^-rich  mixture  could  generate  up  to  10  moles  of  F-atom  (after  due 
allowance  is  made  for  the  enthalpy  change  associated  with  raising  the  gas 
temperature  to  1,  700 “K,  which  corresponds  roughly  to  the  temperature  for 
complete  dissociation  of  F.,  under  thermodynamic  equilibrium  condition  at 
standard  density).  Similarly,  the  heat  of  combustion  from  1  mole  of  Al  in 
a  Cl^-rich  mixture  could  generate  up  to  2.4  moles  of  Cl-atom.  The  corre¬ 
sponding  maximum  F-atom  and  CJ.-atom  yield  from  burning  of  1  mole  (or 
24.  3  grams)  of  Mg  in  an  F  -rich  or  a  Cl  -rich  mixture  is  about  6  and  1.6 


H3 


moles,  respectively.  If  the  l'-atom  or  Cl-alom  so  generated  can  be  utilized 
l'or  direct  pumping  of  such  exothermic:  hydrogen-halide  laser  reactions  as  ’ 


F  4  II  -  HF  i  11  -I  11.7  1.,:?l/moF  (115) 

L, 

F  i  HBr  11F  i  Br  1  45.3  kcal/mole  ( 1 1  tS ) 

F  +  HI  -  HF  -I  I  +  62.  3  kcal/mole  (117) 

Cl  4  HI  -  HC1  +1+32  kcal/mole  (118) 


and  assuming  that  only  12%  of  the  exothermicity  can  be  extracted  as  laser 
energy,  then  combustion  of  1  mole  of  aluminum  could  generate  as  much  as 
38,  55,  and  75  kc.al  of  HF  laser  energy  from  reaction  (115),  (116),  and  (117), 
respectively;  or  4.2  kcal  of  MCI  laser  energy  from  reaction  (1  ifti.  Comparing 
these  with  tlie  ignition  energy  All  ^  4  kcal/mole  indicated  for  Al  in  Table  II, 
it  is  seen  that  a  potential  energy  gain  by  a  factor  of  somewhere  between  2  and 
20  is  quite  possible  through  a  la  i  or -ignition  process  in  a  two-phase  system. 
When  a  conventional  flashlamp  instead  of  a  pulsed  laser  is  used  to  initiate 
the  ignition  process,  the  potential  energy  gain  is  ;vcn  more  significant  since 
the  minimum  energy  required  for  direct  pholodissocialion  of  10  moles  of  F 
is  310  kealories  and  absorption  of  flashlamp  energy  by  solid  particles  is 
likely  to  be.*  more  efficient  than  by  a  dissociating  gas. 

In  order  to  make  use  of  the  heat  of  combustion  from  solid  particles  for 
cd'fji  ieiit  pumping  of  chemical  lasers  in  two-phase  systems,  a  number  of 
physical  and  ihemhal  conditions  must  he  simultaneously  satisfied.  The 


Table  I! 

COMPARISON  OF  IGNITION  ENERGY  AND  HEAT  OF  COMBUSTION  FOR  SOME  SOLID 
PARTICLES.  (Values  taken  from  Refs.  l'.'-  .md  56.  Heat  of  combustion  m 


more  important  of  these  conditions  are:  pH  The  system  must  be  initially  stable, 
or  at  least  quasi- stable  to  the  extent  that  no  significant  spontaneous  conversion 
of  the  reactants  into  products  takes  place  in  the  time  scale  required  to  prepare 
the  two-phase  mixture  for  ignition  within  the  laser  cavity  (a  time  scale  which 
needs  not  be  much  greater  than  the  typical  gas-filling  time  of  a  few  milli¬ 
seconds,  depending  on  the  size  and  geometry  of  the  cavity),  lb)  The  intensity 
of  the  initiating  radiation  source  must  be  suilicicntly  high  so  that  the  solid 
particles  can  be  brought  up  to  ignition  tempera  lure  without  excessive  heat  loss, 
and  that  ignition  of  all  the  particles  will  be.  nearly  simultaneous  in  order  to 
maintain  a  high  degree  of  spatial  homogeneity  for  the:  resulting  laseable  medium, 
(c)  The  burning  rate  of  the  solid  particles  must  be  sujfh  ienily  last  in  compari¬ 
son  with  the  rate  of  de-exritation  of  the  reaction  product  so  that  a  high  popula¬ 
tion  inversion  density  can  be  maintained  for  laser  generation,  id)  The  heat  of 
combustion  generated  in  the  flame  zone  surrounding  the  solid  particles  must 
be  transferred  out,  either  through  radiation  or  through  molecular  transport, 
into  the  ii? t  erpartic)  c  space  and  redistributed  among  Hie  gas-phase  reactants 
in  a  time  scale  shorter  than,  or  at  most  comparable  to,  the  characteristic 
time  for  do -activation  of  the  resultant  gas-phase  reaction  products,  (e)  The 
processes  of  combustion  energy  redistribution  and  subsequent  gas-phase 
reactions  must  favor  the  production  of  the  upper  laser  state  over  that  of  the 


lowei  states. 


It  may  be  noted  that  condition  (a)  is  actually  a  prerequisite  for  any  true 
(i.e.,  premixed)  pulsed  chemical  laser  operation,  and  the  satisfaction  of  which 


is  just  a  matter  of  finding  compatible  mixtures.  Satisfaction  of  condition  (e) 
is  also  somewhat  automatic  if  the  gas-phase  reaction  chosen  were  one  of  the 
principal  pumping  reactions  in  single-phase  chemical  lasers  like  reactions 
(115)  through  (118),  and  that  redistribution  of  the  heat  of  combustion  favors 
the  generation  of  active  atoms  which  participate  in  such  principal  pumping 
reactions.  Condition  (Is)  is  common  to  all  initiation  process,  and  the  satisfac¬ 
tion  of  which  depends  only  on  the  availability  of  an  initiation  source  of  suffi¬ 
ciently  high  intensity. 

Even  though  much  of  the  basic  chemical  kinetics  information  required 
for  quantitative  analyses  of  processes  related  to  the  satisfaction  of  conditions 
(c)  and  (d)  is  not  yet  available,  ordcr-of- magnitude  arguments  based  on 
extrapolation  of  existing  knowledge  about  metal  powder  combustion  in  air  or 
in  other  oxygen-laden  gas  mixtures  does  not  indicate  that  these  conditions 
cannot  be  satisfied.  In  fact,  if  the  unknown  rate  constants  for  exothermic  gas- 
phase  reactions  of  the  type 

A1  +  F  A1F  +  F  i  117.5  kcal/mole  (119) 

Li  1 1  I  Cl2  -  HC1  1  Li  Cl  +  108.  0  kcal/mole  (120) 

were  anywhere  comparable  to  those  for  the  halogen-hydride  reactions  (1) 
through  (4)  cited  in  the  introductory  section,  one  could  even  become  optimistic 
about  the  prospect  of  efficient  two-phase  chemical  laser  operation.  We  will 
demonstrate  this  point  through  the  following  discussion  of  various  critical 
problems  related  to  such  prospects. 


I 


HI 


i 


(?,)  Chemical  an'' 


liovmal  Stability  ot  Two  -  pha  s  c  Systems.  If  the 


mixing  of  finely-dispersed  metal  powders,  such  as  A1  and  Mg,  with  active  gas 
mixtures,  such  as  F  -II  Cl  -HI,  Cl  -HBr,  F  0-11  ,,  etc,,  were  strongly 
hypcrgolic  (i.  e.  ,  spontaneous  ignition  temperature  equal  to  or  below  room 
temperature,  with  very  rapid  combustion  rate),  no  ignition  system  will  be 
necessary.  For  such  hypergolrc  mixtures,  the  only  possible  mode  of  pulsed 
laser  operation  will  be  that  of  an  intermittent -flow  (>,e.  ,  blow-down  type) 
mixing  laser  in  which  the  reactive  components  are  instantly  introduced  into 
tne  laser  cavity  in  much  the  same  way  r  a  reactants  are  introduced  in  any 
single-phase  mixing  laser  system.  Tie  efficiency  oJ  such  Lwo-pliase  mixing 
laser  system  will  be  governed  by  a  competition  among  the  rates  of  reactants 
mixing,  combustion,  halogen- atom  generation,  acti.v.  molecule  formation 
and  colli sional  dc -excitation  of  the  active  molecules. 

If  the  two-phase  mixture  won  stable  or  only  weakly  hypcrgolic  (in 
the  sense  that  the  rates  of  reaction  at  room  temperature  remains  so  low  that 
no  significant  conversion  of  reactants  into  products  takes  place  in  the  time 
scale  required  to  prepare  the  two-phase  mixture  and  to  send  it  through  the 
laser  cavity),  genuine  pulsed  operation  with  premixed  reactants  would  be 
possible.  Economy  of  the  initiation  process  would  then  be  governed  by  the 
specific  heat  of  the  solid  particles  and  the  ignition  temperature  of  such  particles 
in  the  active  gas  mixture.  The  efficiency  for  subsequent: conversion  of  chemical 
energy  into  laser  energy  is  then  governed,  by  a  competition  among  the  rates 


8K 


of  combustion,  active-atom  generation  and  diffusion,  excited  molecule  forma¬ 


tion,  and  collisional  de-excitation  of  the  active  molecules. 

(3)  Kate  of  Heating  of  Solid  Particles.  To  minimize  the 
energy  required  for  bringing  the  surface  temperature  of  the  solid  particles  up 
to  ignition  temperature  without  grossly  heating  the  gaseous  components  in  a 
two-phase  system,  it  would  be  desirable  to  heat  up  the  solid  particles  as 
rapidly  as  possible,  i.  e.  ,  in  a  time  scale  that  is  short  in  comparison  with  the 
characteristic  time  t  ^  for  diffusion  of  heat  in  the  gas  phase  over  the 
averaged  inter-particle  distance  i  (see  subsection  c  above). 


-> 


The  thermal  diffusivity  D  is  generally  a  function  of  gas  density,  temperature 

57 

and  chemical  composition.  However,  for  the  purpose  of  the  present  discus¬ 
sion,  it  suffices  to  consider  D  to  be  identical  to  the  averaged  molecular  dif- 
fusivity  with  a  typical  value  D  -  0.  25  em^/sec  for  all  gas  mixtures  of 
interest  at  standard  density  pn  (defined  here  as  the  donsitv  at  which  the  total 

V' 

19  3 

number  density  of  molecules  in  the  gas  phase  is  equal  to  2.  5  X  10  /cm  ). 
Thus,  at  any  given  gas  density  p,  the  characteristic  interparticle  diffusion 
time  is  given  by  (when  1  is  measured  in  microns), 

i  ^  p  -  8  2  n 

t  —4X10  t  —  second  (122) 

u  uo  po  p0 

It  is  seen  that  for  a  typical  inter-partic.le  distance  of  50  microns,  the 
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characteristic  diffusion  time  is  about  100  ;isec  at  standard  density. 


The  rate  of  heating  of  very  small  particles  by  a  suddenly  turned  on 
source  of  electromagnetic  radiation  was  discussed  earlier  in  subsection  b. 
and  the  calculated  rate  of  temperature  rise  per  unit  radiation  intensity  was 
plotted  in  Fig.  7.  As  noted  earlier,  for  the  healing  ot  aluminum  particles,  in 

-  3 

order  to  obtain  a  temperature  rise  of  1,  0()0'“K  in  a  time  scale  of  10  sec.  , 

6  ?. 

an  intensity  of  10  watls/cnT  would  be  noe»ied  at  1  laser  wavelength,  or 
8  2 

10  watts/cm  would  be  needed  lor  10.6  radiation.  For  the  heating  of 

gxaphite  particles,  the  intensity  requirements  lor  providing  the  same  rate  cf 

3  j 

temperature  rise  is  reduced  to  2.  2  x  10  watts /em*”  at  =  1  u  and 


,5 


i:  Y  y"  tili'jlj  S.ppc-c.1  i* 


to  be  well  within  current  puised  laser  technology.  For  example,  a  c.om- 

merically  available  glass  laser  with  a  9 -joule  energy  output,  30-nsec  pulse 

width,  and  a  beam  diameter  of  12-mm  would  provide  a  p(>ak  intensity  of  about 

8  2 

3  X  10  watts/cm  at  1,06-u  wavelength.  The  same  range  of  intensity  at 
10.6-u  wavelength  can  also  be  reached  by  mild  focussing  of  the  output  beam 
from  Transverse  Electrically-excited  Atmospheric-pressure  (TEA)  CO 

L-. 

lasers  of  performance  characteristics  comparable  to  those  reported  in  the 

..  58,59 

recent  literature. 


K.  g.  ,  American  Optical  Corp.  Model.  No.  A0-30P  as  listed  in  the  Laser 
Focus  1971  buyers’  Guide,  pp.  120-121. 
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When  a  flashlamp  instead  of  a  pulsed  laser  is  used  as  the  initiation 

source,  the  maximum  intensity  available  for  solid  particle  heating  will  be 

limited  by  the  apparent  brightness  temperature  of  the  flashlamp.  Thus,  for 

a  flashlamp  operating  at  a  peak  brightness  temperature  of  15,  000  °K,  the 

5  Z 

maximum  intensity  will  be  limited  to  about  3X10  watts /cm  .  However, 

since  the  light  output  from  such  a  flashlamp  would  have  a  broad  spectral 

distribution  which  peaks  around  3,  000  A,  it  appears  that  such  an  intensity 

would  be  quite  adequate  for  heating  up  highly  conductive  metal  powders, 

such  as  aluminum  and  magnesium,  to  their  respective  ingition  temperature 

-  5  * 

in  a  time  scale  of  the  order  of  10  sec. 

(4)  Burning  Rate  of  Metal  Powders.  The  rates  of  combustion  of 
metal  powders,  such  as  Al,  Mg,  B,  Ti,  Zr,  etc.,  in  gases  such  as  F  Cl 

t-i  C* 

NF ^ ,  F^O,  CO,.,,  etc.  ,  are  generally  not  known.  Hence,  one  can  only  infer 
such  rate  from  experimental  data  published  in  the  literature  for  burning  of 
metal  powder  in  air  and  in  oxygen- laden  gas  mixtures. 

According  to  the  recent  paper  of  Wilson  and  Williams,  ^  the  charactc 
istic:  time  for  combustion  of  aluminum  particles  in  oxygen-laden  gases 

from  all  published  experiments  can  be  reproJe*  ted  by  the  formula 


It  may  be  noted  that  in  order  to  calculate  the  heating  rate  by  a  flashlamp 
properly,  equations  should  be  suitably  modified  to  take  into  account  the  broad 
spectral  distribution  of  the  flashlamp,  as  well  as  any  resonance  absorption 
effect  near  X  =  2iTa. 


VI 


(123) 


r>  Pp(o,) 


where  a  is  the  initial  radius  ol  the  aluniinuin  particle,  p(G  )  is  the  particle 
pressure  of  expressed  in  units  of  the  standard  atmosphere,  and  f?  is  a 
constant  which  depends  only  on  the  dilution  ratio  of  the  oxygen-laden  gas 
mixture  for  any  given  diluent.  For  combustion  of  aluminum  particles  of 

initial  radius  ranging  between  H)  and  300  micron  in  0.,/Ar  mixtures,  the 

-  3  .1  -  4  2. 

value  of  3  was  found  to  be  4  xlO  cm  "/sec  tor  pure  O  and  9  X  10  cm  / 

Ci 

sec  for  a  10°, ’0  0^-90%  Ar  mixture  of  total  pressure  ranging  between  1  and 

5  atmospheres.  The  characteristic  time  t  is  defined  hero  as  the  time  it 

h 

takes  the  combustion  process  to  reduce  (he  initial  radius  of  the  aluminum 
particle  by  a  factor  of  2  (and  hence,  its  volume  by  a  factor  of  8).  Thus, 
according  to  these  experimental  data,  the  burning  time  for  aluminum 
particles  of  1  0  micron  initial  radius  in  pure  oxygen  at.  standard  atmospheric 
pressure  is  about  Z.  i  XlO  ^  sec. 

bO  l 

As  noted  by  Wilson  and  Williams,  the  (X  dependence  of  the 
characteristic  burning  time  T  as  indicated  by  Eq.  (133)  is  consistent  with 
what  one  would  predict  from  the  classical  quasi- steady  vapor-phase  com¬ 
bustion  theory  (sec  Fig.  11)  in  which  the  diffusion  flame  geometry  is  assumed 
self- similar  and  the  rate  of  release  of  combustion  energy  from  the  reaction 


zone  is  assumed  to  be  governed  by  local  equilibrium  chemistry.  Since  the 


chemical  kinetic  rates  never  entered  such  a  theoretical  model  of  solid  particle 
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combustion,  there  is  no  size  limit  on  the  applicability  of  tile  theory  as  long 
as  the  molecular  mean-frec-petli  remains  sufficiently'  small  in  comparison 
with  the  particle  diameter  so  that  the  continuum  treatment  was  justified. 
However,  from  elementary  considerations,  it  would  be  unrealistic  to  expect 
that  the  empirical  formula  (133)  could  be  extrapolated  indefinitely  for  pre¬ 
diction  of  the  burning  time  for  very  small  partn  los  (i.c.  ,  a  ^  10  u  )  flue 

to  the  eventual  onset  of  finite-rate  chemistry  and  lice  molecule  flow.  (Note 
that  the  viscosity  mean- free-pat h  for  air  at  standard  density  is  about.  0.  On  u,  ) 
At  this  juncture,  if  one  assumes  that  the  rates  of  combustion  of  metal 
powders,  such  as  Al,  Mg,  1>,  Ti,  X  r,  etc.,  in  gases  suciias  K  ,,  Cl,,  1-  ,0, 
etc.,  are  at  least  as  fast,  as  that  found  for  combustion  of  aluminum  in  O., 
and  that  the  partjcle-sixe-dcpenden<_e  of  the  charade  ristic  burning,  time  as 
indicated  by  Eq.  (1.13)  can  lie  extrapolated  dow  n  to  a()  of  the  order  of  1  c: , 
then  metal  powder  burning  time  of  the  order  ot  10  ^  sec  appears  reachable 
in  two-phase  ignition  systems  operating  at  standard  density. 

(5)  Redistribution  of  Combustion  Energy  tor  Laser  Humping,  for 
efficient  conversion  of  combustion  energy  into  laser  energy  in  two-phase  sys¬ 
tems,  it  is  essential  that  the  heat  of  combustion  generated  m  the  flame  zone 
surrounding  the  solid  particles  be  transferred  out  into  the  interpu  rtiele  space 
and  redistributed  among  the  gas-phase  reactants  in  a  time  scale  which  is 
shorter  than,  or  at  most  comparable  to,  the  characterisin'  time  for  colli¬ 
sions!  de -activation  of  the  lasing  molecules.  In  chemical  systems  which 
make  use  of  dissociated  atoms  for  pumping  of  the  upper  laser  states  through 
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reactions  of  the  type  represented  by  Reactions  (115)  through  (118),  re¬ 
distribution  of  the  combustion  energy  must  result  in  a  preferential  dissocia¬ 


tion  of  the  halogen-rich  gases  suc  h  as  F  ,,  Cl )(  NF  F  O,  etc.  over  that 

c  J 

of  the  hydrogen- containing  molecules  sucli  as  Ii^,  IIBr,  HI,  etc.  Alter- 
natively,  it  redistribution  of  the  combustion  energy  were  to  result  in  a  pref¬ 
erential  dissociation  of  the  hydride  molecules,  the  JI  atoms  so  generated  may 
be  utilized  for  pumping  of  hydrogen -halogen  lasers  through  "hot  reactions" 
of  the  type 

II  t  F -  ILF  4  F  4  101.5  kcal/mole  (124) 

H  i  Cl,  -  1  TCI  4  Cl  -  45.  1  kcal/mole  (125) 

The  worst  situation  one  can  imagine  is  that  re-distribution  oi  the  combustion 
energy  leads  to  simultaneous  decomposition  of  all  gas-phase  reactants  at 
about,  the  same  rate,  but  such  a  situation  is  probably  not  very  likely. 

Redistribution  of  the  combustion  energy  from  the  hot  flame  zone  cut 
to  the  cold  gas  occupying  the  bulk  of  the  interpariiclc  space  will,  proceed 
radiatively  as  well  as  through  molecular  transport.  If  the  flame  zone  were 
strongly  chemiluminescent  and  the  omitted  light  acre  to  have  a  spectral  range 
which  overlaps  that  for  photodissocialn  of  the  halogen-rich  molecules,  then 
direct  conversion  of  a  significant  fraction  of  the  combustion  energy  into  dis¬ 
sociated  halogen  atoms  could  take  place  instantaneously  throughout  the  inter¬ 
particle  space.  In  other  words,  the  burning  metal  powder  could  act  like  an 
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"internal  ilashlamp"  for  generation  of  active  atoms  in  much  the  same  way 
that  active  atoms  are  generated  by  external  light  sources  in  the  conventional 
flashlamp  initiation  process  for  single -phase  chemical  lasers. 

Visible  light  generated  troni  burning  o)  line  magnesium  wires  or 

powder  in  oxygen  (or  air)  has  lung  been  utilized  as  a  convenient  light  source 

in  flash  photography.  Likewise,  one  may  expect  that  tin  blue  and  near- 

ultraviolet  end  ot  the  emission  spectrum  iron  the  flame  zone  associated  with 

burning  of  aluminum  and  magnesium  powders  in  halogen-rich  gases  can  also 

be  utilized  for  pholodis  socia  t  ion  oi  ]■' ,  and  t.l  ,,  which  are  known  to  have 

strong  dissociation  eontimia  in  the  sped  rut  range  d,  Out)  x  ’■  1,  Oh'O  A. 

(  Photodissoeiation  of  11  from  its  ground  electronic  stale  i  i  not  very  likely 

since'  tile  dissociation  continuu  In  in  fit.  u(r  u  ’  t  ra  \  bbet .  i  According  to 

5  9 

Pearsc  and  Ciaydon,  All'  and  AHd  arc  known  to  have  strong  electronic 
hands  emitting  in  this  spectral  range  of  interest,  hi  addition,  eontimia  emis¬ 
sion  resulting  from  molecular  complex  formation  such  as  fAlF  ,)  ,  (AK.'h  )  , 

3  n  d  n 

or  from  condensation  of  such  reaction  products  into  very  small  refractory 
solid  particles  near  the  high-tempo  rain  re  llunie  zone  may  well  extend  into 
this  near  ultraviolet  region.  However,  one  must  say  at  this  time  thai  we  do 
not  yet  know  enough  about  the  probability  of  excited  electronic  state  forma¬ 
tion  within  the  flame  zone  of  burning  metal  powders,  nor  do  w  e  know  the 
oscillator  strengths  associated  with  such  excited  electronic  states  of  reaction 
products  so  formed,  to  make  any  quantitative  prediction  about  the  efficiency 
of  such  photo  processes.  For  this  reason,  we  suggest  that  some  shock-tube 
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experiments  be  carried  out  to  collect  quantitative  information  on  chem¬ 
iluminescence  from  burning  metal  powders  that  would  be  useful  lor  further 
analysis  of  this  particular  problem. 

In  contrast  to  the  photodissociation  processes,  which  can  be  con¬ 
sidered  instantaneous  for  all  practical  purposes,  redistribution  of  combus¬ 
tion  energy  into  the  bulk  of  the  interparticle  space  through  molecular 
diffusion  and/or  hoat  conduction  would  be  much  slower.  Even  though  the 
detailed  kinetic  steps  related  to  the  burning  of  metal  powders  are  not  known 
(not  even  for  burning  of  aluminum  powder  in  oxygen,  according  to  experts  in 
the  field  of  combustion  theory*1*  hone  may  postulate  that  gas-phase  combus¬ 
tion  of  metal  vapors,  such  as  A1  and  Mg,  in  diatomic  gases  such  as  and 

i — ' 

Cl  ,  is  likely  to  proceed  according  to  sequences  of  the  type 


A1 

-L 

r* 

‘  2 

F 

4 

117.5  kcal/mole 

A1F 

+ 

■l  2 

F 

4 

8  7.  1  kcal/mole 

A1F2 

1 

F2 

F 

4 

1  06.  0  kcal/mole 

A1 

+ 

Ci2  - 

A1C1 

3  Cl  4  5  3.  1  kcal/mole 

A1C1 

+ 

C12  - 

AlCi 

f  Cl  i  35,  3  kcal/mole 

A1C1 2 

+ 

CI2  - 

A1C1 

3 

+  Cl  3  35.  5  kcal/mole 

(126) 


(127) 
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Out  of  the  310.6  kcal/mole  generated  from  sequence  (136),  or  the  133.  9 
kcal/mole  generated  from  sequence  (127),  only  70.4  kcal/mole  would  be 
needed  for  vaporization  of  the  aluminum  solid.  The  remaining  240.  2  kcal/ 
mole  from  sequence  (126)  or  53.5  kcal/mole  from  sequence  (127),  would 
either  be  radiated  out  into  the  inter-particle  space  or  deposited  in  the  reac¬ 
tion  products  A1F,  A1F  ?,  AIF^,  F,  or  A1C1,  A1C1  A1C1  ,  Cl,  in  the  form 
of  vibrational  and  translational /rotational,  energies.  In  the  case  of  sequence 
(126),  outward  diffusion  of  the  energetic  reaction  products  into  the  inter- 
particle  space  may  well  lead  to  further  dissociation  of  the  ambient  F  through 
reactions  of  the  type 


AlKlv)  -t  F,  ^ 

Li 

ai  K(o}  3  r 

3  ]• 

(128) 

A1F  ,(v)  1  F  ,  - 

All’,(0)  l  F 

3  1  ‘ 

(129) 

A1F  (v)  3  F  - 

A1F  (0)  3  V 

3  F 

(130) 

F  •«  F2  - 

F  3  F 

3  ]•• 

(131 ) 

For  sequence  (127!,  however,  the  residual  eueigy  from  none  of  the  three 
steps  appears  high  enough  for  further  dissociation  of  Ci  but.  by  the  end  of 
step  3,  three  moles  of  atomic  Cl  would  already  have  been  produced  from 
burning  of  one  mole  of  aluminum. 

(6)  Matching  of  Molecular  Diffusion  Time  and  Collisional  Deexcita¬ 
tion  time.  If  chemiluminescence  were  weak  and  (lie  photodissociation 
processes  were  completely  ineffective  in  redistributing  the  combustion 
energy,  one  must  then  rely  on  diffusion  of  the  active  alums  and/or  the 
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energetic  molecules  for  such  energy  redistribution. 

Again,  going  back  to  the  example  of  the  hydrogen-halide  systems  cited 
in  the  preceding  section,  one  may  readily  note  that,  concurrent  with  the  dif¬ 
fusion  of  the  F  or  Cl  atoms  into  the  interparticle  space,  laser  pumping  reac¬ 
tions,  such  as  Reactions  (115)  through  (118)  will  take  place.  The  excited 
laser  molecules,  such  as  lll'(v)  or  )TCl(v),  which  are  formed  earlier  (i.e,, 
closer  to  the  flame  zone)  will  also  be  .subjected  to  more  de -excitation  colli¬ 
sions  as  the  outward  diffusion  of  active  atoms  and  energetic  molecules 
progresses.  Thus,  in  order  to  maintain  a  positive  gain  over  the  bulk  of  the 
into r -particle  volume,  it  is  essential  that  the  characteristic  diffusion  time 
~D  l3C  shorter  than,  or  at  most  comparable  to,  the  characteristic 

collisional  deexcitation  time 


(132) 


where  k_  denotes  some  averaged  rate  constant  for  collisional  de- excitation 


of  the  upper  k 


•  iw,  lu  u  y  c*.  v.  y  \j^—  |  liiwumui  v; 


1  ecu le  over  lIic  mlerpa  r  i  icli; 


space. 


in  denotes  the  corresponding  averaged  number  density  of  type-j  molecules, 
and  the  summation  is  to  be  carried  over  all  types  of  molecules  present  in 


the  gas  mixture  filling  the  in terparticle  space.  By  writing 


n.  =  2.  5  X  1 0 1  'J  I  M.  ]  —  cm 

1  P0 


(133) 
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whore  [M  ]  denotes  some  averaged  mole  fraction  of  type-i  molecules 
(appropriately  normalized  with  respect  to  the  initial  molar  concentratioxx  of 
the  total  gas-phase  mixture  before  combustion  and  dissociation  took  place), 
and  p  / p  is  the  actual  gas  density  divided  by  the  standard  gas  density  as 
defined  earlier.  The  characteristic  collisional  de-excitation  time  may  be 
expressed  in  the  form 

1 


I  •>  c  P  NT  .w  1  I 


I 


Z.  5  X  10 


0 


2  \.j 


) 


second 


(134) 


if  k  .  is  given  in  units  of  cm  /molecular- sec. 
cj 

By  equating  Eq.  (122)  and  Eq.  (134),  one  therefore  obtains  the  follow¬ 
ing  expression  for  matching  the  characteristic:  diffusion  time  T  with  the 
characteristic  collisional  de-excitation  time  t 


■ 10-6  |2kcj|Mj' 


-i/. 


(135) 


where  t  is  the  averaged  interparticle  distance  measured  in  microns.  The 

collisional  de-excitation  rate  constant  k  .  depends,  of  course,  strongly  on 

v  J 

the  type  of  laser  molecule  and  the  particular  upper  laser  state  under  consid¬ 
eration.  It  may  also  depend  strongly  on  the  chemical  type  of  the  collision 
partner  Ivh  ,  and  on  the  translational/ rotational  temperature.  For  example, 

the  value  of  k  .  for  collisional  de-excitation  of  the  v  =  1  level  of  HF  at 
cj 

- 17  3 

T  -  500°K  is  about  8X10  cm  /molecule-sec  when  the  collision  partner 
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-12  3 

is  Ar  or  F  ;  and  about  2.  5  XlO  cm  /molecule- sec  when  the  collision 

partner  is  another  HF,  At  a  temperature  of  1,  500DK,  the  corresponding 

- 15  3 

value  of  k  .  becomes  7  X  10  cm  /molecule- sec  for  collision  with  Ar 
cj 

-12  3 

or  F  and  1.7x10  '  cm'  /molecule-sec.  for  collision  with  another  HF, 

respectively.  On  the  other  hand,  the  value  of  k  .  for  collisional  de- 

cj 

- 14  3 

excitation  of  the  v  -  1  level  of  HC1  at  T  -  500lK  is  about  1.  7  XlO  cm  / 


molccule-sec,  whether  the  collision  par  tner  were  Cl ^  ,  H  H,  or  another 

HC1  molecule  at  ground  vibrational  level.  At  a  temperature  of  1,  500°K, 

- 1  3  3 

the  value  of  k  .  becomes  4  XlO  cm  /molccule-sec  for  the  same  set  of 

cj 

collision  partners. 


-  1  2  3 

By  taking  typical  value  £k  .  [M.]  -  10  cm  /molecule -  sec  fo  r 

CJ  } 

- 1  4  3 

collisional  de-excitation  of  HF  vibrations,  and  Id  cm  /molecule-sec  for 


collisional  de- excitation  of  HC1  vibrations,  respectively,  Eq.  (1 35)  indicates 
that  matching  of  t  and  "  can  be  achieved  for  J  p/p  2  1  micron  in  an 

i—  ’  W  V 

HF  system,  and  for  £  p/Pq  -  10  micron  in  an  IICl  system.  Thus,  at  standard 

density,  the  characteristic  diffusion  time  t.  can  be  made  shorter  than  the 

JJ 

characteristic  collisional  de-excitation  time  t  if  the  averaged  inter- 

v. 

particle  distance  f  is  kept  shorter  than  about  1  micron  in  an  IIF  system, 


and  10  microns  in  an  HC1  system.  At  one- tenth  the  standard  density  (i.  e.  , 


P/Pq  ::  0.  1,  the  corresponding  value  icr  l  should  he  kept  shorter  than  10 
micron  in  an  HF  system,  and  100  microns  in  an  HC1  system,  in  order  to  make 
the  diffusion  time  shorter  than  the  collisional  de-excitation  time  of  4X10  ^  sec 
for  HF,  and  4  XlO  ^  see  for  HC1,  respectively . 
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(7)  Simultaneous  Matching  of  All  Necessary  Conditions,  After 
matching  the  molecular  diffusion  time  with  the  coliisional  de¬ 

excitation  time  T  in  accordance  with  Eq.  (135),  one  must  now  go  back  and 
re-examine  all  other  constraints  that  may  ha.ve  been  imposed,  either  directly 
or  indirectly,  on  the  product  f  p/p  through  the  requirement  for  matching 
of  other  characteristic  time  scales  discussed  earlier,  and  then  see  if  they 
are  all  compatible. 

The  averaged  interparticle  distance  9  is  related  to  the  number  of 

solid  particles  per  unit  volume,  n  ,  such  that  i  =  n  ^  ^  (since  the  averaged 

volume  occupied  by  each  particle  is  n  ^  ).  The  initial  volume  fraction  of 

solid  narticles  <t>  is  niven  bv  n  times  the  initial  volume  ner  nartjcle 
*  ~  '  p  '  p  * 

4  tr  / 3 ,  so  that 


*T  (t1 


(136) 


The  volume  fraction  of  solid  particles  0^  to  be  used  in  a  two-phase  system, 
in  turn,  is  related  to  the  mole  ratio  between  the  solid-phase  particles  and 
the  gas-phase  reactants,  such  that 


r  _  No,  of  moles  of  solid  particles  per  unit  volume 
p  No.  of  moles  of  gas-phase  mixture  per  unit  volume 


molal  volume  of  gas-phase  mixture  at  standard  density 
(P/Pq)  molal  volume  of  solid 


2.  24  X  1 0  0 

T7  V  V"1 


(137) 


10Z 


By  combining  (136)  and  (137),  one  obtains 


=2.  2X  10"Z  (V  X  )1/3(p/p  )1/3 
i  I1  P  0 


(138) 


The  rnolal  volume  V  for  most  metals  of  interest'  to  the  present  problem 

3  3 

is  of  the  order  of  10  cm  (e.  g.  ,  V  -  10.  0  cm  f or  Al,  1 4,  0  for  Mg,  10.6 
for  Ti,  14.  3  for  Zr,  4.  3  for  13,  ek  .  ),  while  the  desired  mole  ratio  X 

p 

between  solid  and  gas  is  expected  f<>  be  of  the  order  0.  1.  Therefore,  one 
may  let  (V^X^)  ^  —  1  for  all  practical  purposes.  The  required  ratio  be¬ 
tween  (.he  initial  particle  radius  and  the  averaged  interparticle  distance  is 
the r e f o r o ,  app r oxinia t e ly 


=  2.  2  >'-  JO  (p/p  )^i 


(139) 


Thus,  if  one  denotes  the  value  of  /  p/p y  which  satisfied  the  matching  condi- 

* 

tion  (135)  by  the  symbol  (t  p/p  )  ,  Kq.  (139)  requires  that  the  initial  radius 
of  the  particles  must  satisfy  the  following  inequality  in  order  to  keep 
T  £  T 

D  C 


CL  £  2,  2  X  1 0 


-  2  ('P/P0’ 


(140) 


In  the  preceding  section,  we  estimated  that  (.(  p/p  )  =1  micron  for  the  HF 

system,  and  (i  p/p  )  ~  10  iTiicron,  Accordingly,  Eq.  (140)  requires  that 

-  2 

the  particle  radius  be  kept  smaller  than  Z.  2.  XlO  and  0.  22  micron  for 


normal  density  operation  in  HF  and  HC1  systems,  re  spec. tively.  Foi  opera¬ 
tion  at  P/Pq  -  0.  1,  the  initial  particle  radius  requirement  becomes 
aQ  ^  0.  102  and  1.  02  micron  for  HF  and  1IC1  systems,  respectively. 

If  the  particle  burning  time  formula  (123)  can  be  extrapolated  to 

burning  of  aluminum  particles  down  to  the  sub-micron  region  in  F^  and  Cl  , 

-  3  2 

then  using  the  value  of  3  -  4  XlU  cm  /sec  previously  obtained  for  Al-O 

-9  -7 

combustion,  one  obtains  a  burning  time  i'  -  1.2X10  and  1.  2  X  10  sec 

-  2 

for  burning  of  2,  2  XlO  and  0.  22  micron  particles,  respectively,  at  normal 

density.  For  burning  at  0.  1  normal  density,  the  burning  time  so  obtained 

-7  -  5 

becomes  2.  6  X10  and  2.  6  XlO  sec  for  burning  of  0.  102  and  1.  02  micron 

radius  aluminum  particles,  respectively.  Since  the  collisional  de-excitation 

-  7  -  5 

time  at  P/Pq  =  0.  1  has  been  estimated  to  be  4  XlO  sec  and  4  XlO  sec 
for  HF  and  HC1  respectively,  these  estimated  burning  times  appear  quite 
satisfactory.  However,  as  we  discussed  earlier,  we  do  not  really  know  if 
the  burning  time  formula  (123)  can  be  extrapolated  to  burning  of  aluminum 
powder  in  halogen  gases,  nor  do  we  know  how  far  it  can  be  extrapolated  in 
the  direction  of  decreasing  particle  radius  beyond  those  covered  by  previous 
experiments. 

From  preceding  discussion  of  the  rate  of  heating  of  solid  particles 
in  subsection  (3),  we  concluded  that  heating  time  scale  of  the  order  of  10  ^  sec 

appeared  quite  feasible,  but  further  compression  of  the  heating  time  scale 

-  5 

to  values  much  shorter  than  10  sec  would  cause  undue  strain  on  the 
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requirement  of  radiation  source  intensity  (unless  the  absorbing  particles 
were  made  of  much  poorer  conductors  than  aluminum).  This  additional 
constraint  accordingly  further  limits  the  feasible  range  of  two-phase  pulsed 
hydrogen-halide  lasers  to  initial  gas-phase  density  somewhat  below  the 
standard  atmospheric  density  (say,  p/p  ~  0.1).  However,  comparing  this 
to  the  accessible  density  range  of  existing  single-phase  hydrogen-halide 
systems,  the  present  estimated  feasible,  gas  density  (for  undiluted  gas- 
phase  reactants)  p/p^  ~  0.1  already  represents  a  very  significant  potential 
for  improvement. 


105 


APPENDIX  ) 


RECOMMENDED  RATE  CONSTANTS  FOR  H2-F2 
CHEMICAL  LASER  (From  Ref.  47) 


A  detailed  list  of  reaction  considered  l'or  the  rate  equations  of  the 
first  and  second  vibrational  levels  of  the  HE  molecule  together  with 
selected  rate  constants  is  presented  here  (only  forward  reactions  are 
considered ). 


Excitation  Reactions 


E  +  H. 


4(1) 


HE  (1 )  t  II 


hu-1'"  1°13'1'6/° 


V  4  II, 


4(2) 


111(2)  l  ll 


k  ,  •-  V  -  JO 
4  (2) 


1  PI.  6/ 


H  4  F 


'5(1) 


I1F(1)  l  F 


kSllP’X,“ 


12-2.  4/0 


H  1  F2  I  IE  (2)  I  E 


k5(2)  =  9X  10 


12-2.  4/( 


De  excitation  React  ion  s 

(1)  V-T  transfer 


HE (1  )  4  HE  —  HE (0)  +  HE  k,  =  106-  23  T1’  ??  4  1016t‘ 

ha(l ) 


HF  (1 )  4  H2  I  IF  (0 )  4  112 


HF(2)  4  1TF  1IF(1)  4  HE  k  ... -2k,  , .  , 

6a  (2)  6a(l) 


.1. 43 


106 


6a  (2) 

HF  (2) 

4  HZ 

- - VHF(U  +  HZ 

HF(3) 

+  HF 

k6a(3) 

- 

HF  (2)  4  HF 

HF  (3) 

+  H2 

k6a(3) 

- 

I IF (3)  4  112 

HF  (1 ) 

+  F 

k6b(l) 

- - 

I  IF  (0)  4  F 

k 

HF  ( 2 ) 

+  F 

6b  (2) 

- > 

I  IF  ( 1 )  i  F 

HF  (3) 

+  F 

k6b(3)a 

I IF  (2)  4  F 

HF  (1 ) 

4  FZ 

k 

6c(l) 

- 

HF  (0 )  i  F2 

HF  (2) 

4  FZ 

6c  (2) 

- > 

Hi’(l)  1  FZ 

HF  (3 ) 

4  FZ 

k6c  (3)^ 

HF (2)  4  FZ 

HF  ( 1 ) 

+  H 

k6f(U 

HF (0 )  I  II 

IIF(2) 

-t  H 

k6f(2) 

— - ^ 

HF  ( 1  )  4  I  I 

HF(3) 

+  H 

k6f(3) 

HF (2)  4  11 

k6a(3)  3k6a(l) 


k , 

=  1.5 

xio10 

'6b  (1) 

=  2  x 

k 

6  b  (Z) 

6b  ( 1 ) 

k 

=  3  X 

k 

6b  (3) 

6b  (1) 

k.  ..  . 

=  5.  45  X  10' 

be  (I  .) 

k 

=  2  X 

k  v 

6  c.  (2) 

6c  (1) 

k 

=  3  X 

6c  (3) 

6c  (1 ) 

k _ 

=  T  (1 . 

.5X10 

61(1) 

k 

=  2  x 

k 

61(2) 

61(1) 

k 

=  3  X 

k ,  , 

61(3) 

61(2) 

T 


10-1.  9/9, 


(2) 


v-v  exchange 

HF  ( 1 )  -1-  HF  (1 ) 

k7a(l) 

- H F ( 0 )  4  IIF (2) 

k7.,(1F2,lc8Tl'5' 

IIF(l)  +  HF  (2) 

k7b(l ) 

- IIF (0)  4  IIF(3) 

k7b(i)  ‘  3(5  *  10?  r1'5: 

HF (1 )  -t  HF (3) 

k7c(l  ) 

- HP'  (0)  4  I IF  (4  ) 

k7.-m  =,(2x!i)7  T‘'5; 

k 

HF  (1 )  +  HF  (3 ) 

a  (2 

- — 1 1 1 ' ( 2 )  4  IIP' (2) 

k7a(2)'fc<10  1  ’ 

HF  ( 2 )  4  HF  (0 ) 

k7a(I ) 

- - — a-  HF(1)  4  IIF(l) 

HF(2)  4  HF  (2) 

k7a(2) 

- ] IF ( 1 )  4  HF  (3) 

HF{2)  4  IIP' (3)  - 

k7b(2) 

- - ♦IIP  (1 )  <  HF(4' 

k7bU>  10?  t1'5) 

where  0  -  2.  3  RT/1000  Real  /mole,  R  is  the  universal  gas  constant, 
T  is  the  temperature  in  °K,  and  the  rate  constants  k  are  in  units  of 
cc  ./mo]  e-sec. 
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